Signature Sequence of 
Intersection Curve of Two Quadrics for 
Exact Morphological Classification 



i> ■ 

Changhe Tu a Wenping Wang ' Bernard Mourrain c 
(N ' Jiaye Wang ; 

0^ \ b University of Hong Kong 

C INRIA 

o 
u 

o 

Abstract 

■ We present an efficient method for classifying the morphology of the intersection 

curve of two quadrics (QSIC) in PR 3 , 3D real projective space; here, the term 
morphology is used in a broad sense to mean the shape, topological, and algebraic 
properties of a QSIC, including singularity, reducibility, the number of connected 
components, and the degree of each irreducible component, etc. There are in total 
35 different QSIC morphologies with non-degenerate quadric pencils. For each of 
these 35 QSIC morphologies, through a detailed study of the eigenvalue curve and 
the index function jump we establish a characterizing algebraic condition expressed 
in terms of the Segre characteristics and the signature sequence of a quadric pencil. 
We show how to compute a signature sequence with rational arithmetic so as to 
determine the morphology of the intersection curve of any two given quadrics. Two 
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immediate applications of our results are the robust topological classification of 
QSIC in computing B-rep surface representation in solid modeling and the derivation 
of algebraic conditions for collision detection of quadric primitives. 
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1 Introduction 



Quadric surface, being the simplest curved surfaces, are widely used in com- 
putational science for shape representation. It is therefore often necessary to 
compute the intersection or detect the interference of two quadrics. In com- 
puter graphics and CAD/CAM, the intersection curve of two quadrics needs 
to be found for computing a boundary representation of a 3D shape defined 
by quadrics. In robotics (1271 ) and computational physics (1201 ; |28| ) one often 
needs to perform interference analysis between ellipsoids modeling the shape 
of various objects. There have recently been rising interests in computing the 
arrangements of quadric surfaces in computational geometry (0; i|), a field 
traditionally focused on linear primitives. 



The intersection curve of two quadric surfaces will be abbreviated as QSIC. 
Exact determination of the morphology of a QSIC is critical to the robust 
computation of its parametric description. We study the problem of classifying 
the morphology of a QSIC in PR 3 (3D real projective space); here, we use 
the term morphology in a broad sense to mean the shape, topological, and 
algebraic properties of a QSIC, including singularity, the number of irreducible 
or connected components, and the de gree of each irreducible component, etc. 
There are many types of QSIC in PR 3 (1321 ) . A nonsingular QSIC can have zero, 
one, or two components. When a QSIC is singular, it can be either irreducible 
or reducible. A singular but irreducible QSIC may have three different types of 
singular points, i.e., acnode, cusp, and crunode, while a reducible QSIC may 
be planar or nonplanar. A planar QSIC consists of only lines or conies, which 
are planar curves, while a reducible but non-planar QSIC always consists of a 
real line and a real space cubic curve. Among planar QSICs, further distinction 
can be made according to how many of the linear or conic components are 
imaginary, i.e., not present in the real projective space. 

There are mainly three basic problems in studying the morphology of a QSIC: 

1) Enumeration: listing all possible morphologically different types of QSICs; 

2) Classification: determining the morphology of the QSIC of two given quadrics; 

3) Representation: determining the transformation which brings a given prob- 
lem QSIC into a canonical representative of its class. We emphasize on the 
second problem of classification, which is an algorithmic issue, while also hav- 
ing the first problem solved as a by-product of our results. Specifically, we 
enumerate all 35 different morphologies of QSIC, and characterize each of 
these morphologies using a signature sequence that can exactly be computed 
using rational arithmetic for the purpose of classification. The third problem, 
not handled here, leads to a lengthy case by case study which depends a lot 
on the application behind. 

Consider the intersection curve of two quadrics given by A: X T AX = and B: 
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X T BX = 0, where X = (x, y, z, w) T E PR 3 and A, B are 4x4 real symmetric 
matrices. The characteristic polynomial of A and B is defined as 



f(X)=det(XA-B), (1) 

and /(A) = is called the characteristic equation of A and B. 

The characteristic polynomial /(A) is defined with a projective variable A G 
PM; thus it is either a quartic polynomial or vanishes identically. The latter 
case of /(A) vanishing identically occurs if and only if A and B are two singular 
quadrics sharing a singular point; thus, all the quadrics in the pencil formed 
by A and B are singular. In this case, the pencil of A and B is said to be 
degenerate; otherwise, the pencil is non- degenerate. For example, if A and B 
are two cones with their vertices at the same point, then they form a degenerate 
pencil. When two quadrics form a degenerate pencil, by projecting the two 
quadrics from one of their common singular points to a plane P not passing 
through the center of projection, we reduce the problem of computing the 
QSIC to one of computing the intersection of two conies in the plane P, which 
is a separate and relatively simple problem. For this reason and the sake of 
space, we will not cover this case in the present paper. Hence, we assume 
throughout that /(A) does not vanish identically. 

Our contributions are as follows. We consider a new characterization of the 
QSIC of a pencil, namely the signature sequence, and show how it can be 
computed effectively and efficiently, using only rational arithmetic operations. 
We establish a complete correspondence among the QSIC morphologies, the 
Segre characterization over the real numbers, the Quadric Pair Canonical 



Form (1251 ; 1491 ; I4d ) and the signature sequence, which allows us to derive a direct 
algorithm based on exact arithmetic for the classification of QSIC. Based on 
this correspondence, a simplified analysis of the morphology of different QSIC's 
is described. We obtain a complete table of all the possible morphologies of 
QSIC, with their Segre characterizations, signature sequences and Quadric 
Pair Canonical Forms. These results apply to any quadric pencil whose char- 
acteristic polynomial /(A) does not vanish identically. The case of /(A) = 
leads to the classification of conies in PM 2 , which is not treated here. Tables 
1, 2 and 3 give the complete list of all 35 different types of QSICs in PM 3 with 
non-degenerate quadric pencils. A detailed explanation of these tables is given 
in Section 2.7. 



We stress that this paper is not about afline classification of QSICs, although 
the results of this paper can be used for an implementation of affine classi- 
fication by further considering the intersection of a QSIC with the plane at 
infinity. 

A few words are in order about our approach. Since any pair of quadrics 
can be put in the Quadric Pair Canonical Form, we obtain all possible QSIC 
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morphologies by an exhaustive enumeration of all Quadric Pair Canonical 
Forms, with distinct Jordan chains and sign combinations. For each pair of the 
Quadric Pair Canonical Forms, on one hand, we obtain its index sequence, and 
on the other hand, we determine its corresponding morphology. The derivation 
of the index sequence necessitates the study on eigenvalue curves and index 
jumps at real roots of a characteristics equation, while the determination of the 
QSIC morphology is largely based on case-by-case geometric analysis of two 
quadrics in their Quadric Pair Canonical Forms. Finally, we convert all index 
sequences to their corresponding signature sequences for efficient and exact 
computation. In this way we establish a complete correspondence among the 
QSIC morphologies, Quadric Pair Canonical Forms and signature sequences. 
Overall, the paper is mainly about an algorithm for determining the type of 
an input QSIC. The algorithm itself is very simple, but it is based on a new 
framework of using the signature sequences of different QSICs. Therefore, the 
large portion of the paper is devoted to identifying the signature sequence 
of each of the 35 QSICs, rather than to describing the flow of the simple 
algorithm. 

The remainder of the paper is organized as follows. We discuss related work 
in the rest of this section. Uhlig's method and other preliminaries, including 
a careful study of the eigenvalue curves of a quadric pencil, are introduced 
in Section 2. For an organized presentation, characterizing conditions for dif- 
ferent QSIC morphologies are grouped into three sections: nonsingular QSIC 
(Section 3), singular but non-planar QSIC (Section 4), and planar QSIC (Sec- 
tion 5). In Section 6 we discuss how to use the obtained results for complete 
classification of QSIC morphologies. We conclude the paper in Section 7. 

For a better flow of discussion, in the main body of the paper we will include 
only the proofs of theorems for the first few cases of QSICs, so as to give the 
gist of the techniques employed. The proofs for the rest cases will be given in 
the appendix. 



1.1 Related work 

Literature on quadrics abounds, including both classical results from algebraic 
geometry and modern ones from computer graphics, computer-aided geomet- 
ric design (CAGD) and computational geometry. Classifying the QSIC is a 
classical problem in algebraic geometry, but the solutions found therein are 
given in PC 3 (3D complex projective space), and therefore provide only a 
partial solution to our classification problem posed in PR 3 . Some methods 
for computing the QSIC in the computer graphics and CAGD literature do 
not classify the QSIC morphology completely, while others use a procedu- 
ral approach to computing the QSIC morphology. The procedural approach 
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is usually lengthy, therefore prone to erroneous classification if floating point 
arithmetic is used or leading to exceedingly large integer values or complicated 
algebraic numbers if exact arithmetic is used. 



When the input quadrics are assumed to be the so-called natural quadrics, i.e., 
special quadrics including spheres, circular right cones and cylinders, there are 
several methods that exploit geometric observations to yield robust methods 



for computing the QSIC (I21l;l22l;l30f). However, we shall consider only methods 
for computing the QSIC of two arbitrary quadrics, and focus on how these 
methods classify the QSIC morphology. 

In algebraic geometry the QSIC morphology is classified in PC 3 , the complex 
projective space using the Segre characteristic (j4j). The Segre characteristic is 
defined by the multiplicities of the roots of /(A) = with respect to /(A) as 
well as the sub-determinants of the matrix XA — B. The Segre characteristic 
assumes the complex field, i.e., assuming that the input quadrics are defined 
with complex coefficients, and therefore it does not distinguish whether a root 
of /(A) = is real or imaginary. When applying the Segre characteristic in 
PR 3 , several different types of QSICs in PR 3 may correspond to the same Segre 
characteristic, thus cannot be distinguished. An example is the case where four 
morphologically different types of nonsingular QSICs correspond to the same 
Segre characteristic [1111], meaning that /(A) = has four distinct roots; (see 
cases 1 through 4 in Table 1). 



QSICs in PR 3 , real projective space, are studied comprehensively in (I la; l33l). 
but the algorithmic aspect of classification is not considered. In this paper we 
obtain a complete classification by signature sequences of quadric pencils and 
apply this result to efficient classification of QSICs in PR 3 . 

A well-known method for computing QSIC in 3D real space is proposed by 



Levin ()18l ; Il9l ). based on the observation that there exists a ruled surface in 
the pencil of any two distinct quadrics in PR 3 . Levin's method substitutes a 
parameterization of this ruled quadric to the equation of one of the two input 
quadrics to obtain a parameterization of the QSIC. However, this method does 
not classify the morphology of the QSIC; consequently, it does not produce 
a rational parameterization for a degenerate QSIC, which is known to be a 
rational curve or consist of lower-degree rational components. 

There have been proposed several methods that improve upon Levin's method. 



Sarraga (]29|) refines Levin's method in several aspects but does not attempt to 



completely classify the QSIC. Wilf and Manor (I48l ) combine Levin's method 
with the Segre characteristic to devise a hybrid method, which, however, is 
still not capable of completely classifying the QSIC in PR 3 ; for example, the 
four different types of nonsingular QSICs are not classified in PR 3 . Wang, 
Goldman and Tu (j46l ) show how to classify the QSICs within the framework 
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of Levin's method. DuPont et al (8j) proposed a variant of Levin's method 
in exact arithmetic by selecting a special ruled quadric in the pencil of two 
quadrics, in order to minimize the number of radicals used in representing the 
QSIC; an implementation of this method is described in (1171 ). The methods 
in (I46l ) and (|8|) both adopt a lengthy procedural approach, with no systematic 
approach for a complete classification. 



A different idea of computing the QSIC, again using a procedural approach, 
is to project a QSIC into a planar algebraic curve and analyze this projection 
curve to deduce the properties of the QSIC, including its morphology and 
parameterization. Farouki, Neff and O'Connor (12) project a QSIC to a planar 
quartic curve and factorize this quartic curve to determine the morphology of 
the QSIC. (Note that only degenerate QSICs are considered in (1121) .) Wang, 
Joe and Goldman (|45l ) project a QSIC to a planar cubic curve using a point 
of the QSIC as the center of projection; this cubic curve is then analyzed to 
compute the morphology and parameterization of the QSIC. However, exact 
computation is difficult with this method, since the center of projection is 
computed with Levin's method. 



The work of Ocken et al (26), Dupont et al (jHI^), Tu et al fj36l ) and (37) all use 
simultaneous matrix diagonalization forjxjmputing or classifying the QSIC. 

is not based on any established 



49c |40| ). and the analysis in (1261 ) 



The diagonalization procedure used in 
canonical form, such as the Uhlig form (l25l ; 
is incomplete - it leaves some cases of QSIC morphology missing and some 
other cases classified incorrectly; for example, the case of a QSIC consisting 
of a line and a space cubic curve is missing and the cases where /(A) = has 
exactly two real roots or four real roots are not distinguished. The classification 
by Dupont (0; 0) is based on the Quadric Pair Canonical Form and involves 
criteria such as signature and sign of deflated polynomials at specific roots of 
the characteristic polynomial, leading to a complete procedure to determine 
the type of a QSIC, covering also the case where the characteristic polynomial 
vanishes identically. 



In the above methods some cases of different QSIC morphologies need to 
be distinguished using procedures involving geometric computation, such as 
extracting singular points or intersecting a line with a quadric. Application 
of such procedural methods is not uniform and follows a case by case study, 
which is very specific to the tridimensional problem. 

It is therefore natural to ask if it is possible to determine the morphology of 
a QSIC by checking some simple algebraic conditions, rather than invoking a 
long computational procedure. 

Several arguments are in favor of more algebra. First, a description of the 
configurations of QSIC by algebraic conditions allows us to introduce easily 
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new parameters in our problem. For instance, introducing the time, it has 
direct application in collision detection problems. Secondly, it provides a com- 
putational framework to analyse the space of configurations of QSIC and the 
stratification induced by this classification, that is how the different families 
are related and what happen when we move on the "border" of these fami- 
lies. Moreover, the correlation between the canonical form of pencils and the 
algebraic characterisation can be extended in higher dimension. 

Algebraic conditions have recently been established for QSIC morphology or 
configuration formed by two quadrics in some special cases. The goal here 
is to characterize each possible morphology or configuration using a simple 
algebraic condition, which can be tested or evaluated easily and exactly to 
determine the type of an input morphology or configuration. In related topics, 
a simple condition in terms of the number of positive real roots of the char- 



acteristic equation /(A) is given by Wang et al in (1431 ) for the separation of 



two ellipsoids in 3D afline space. Similar algebraic conditions are obtained by 



Wang and Krasauskas in (1471 ) for characterizing non-degenerate configurations 



formed by two ellipses in 2D afline plane or ellipsoids in 3D afline space. 



As for QSICs, the Quadric Pair Canonical Form form is used in (1361 ) to derive 
simple characterizing algebraic conditions for the four types of non-singular 
QSICs in terms of the number of real roots of the characteristic polynomial; 
however, two of the four types are not distinguished, i.e., they are covered by 
the same condition. This pursuit of algebraic conditions is extended to cover 



all 35 QSICs of non-degenerate pencils in the report (1371 ). which again uses 



the Quadric Pair Canonical Form to derive characterizing conditions in terms 



of signature sequences. The present paper is based on (1371 ) 



Finally, we mention that Chionh, Goldman and Miller (J5|) uses multivariate 
resultants to compute the intersection of three quadrics. 



2 Preliminaries 



2.1 Simplification techniques 



There are two transformations that we will use frequently to simplify the 
analysis of a QSIC. Based on Quadric Pair Canonical Form results (1251 ; |49 ; 



40) (see also Section 2.3), we sometimes apply a projective transformation 
to both A and B to get a pair of quadrics A' : X T (Q T AQ)X = and B' : 
X T (Q T BQ)X = in simpler forms. The transformed quadrics A' and B' are 
projectively equivalent to A and B, therefore have the same QSIC morphology 
in PR 3 and the same characteristic equation as the pair A and B. 
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Table 1 

Classification of nonplanar QSIC in PR 3 



[Segre] r 
r = the # 
of real roots 


Index 
Sequence 


Signature Sequence 


Illus- 
tration 


Representative 
Quadric Pair 


[link 


JJ 

(1|2|1|2|3> 


(1,(1,2),2,(1,2),1,(1,2),2,(2,1),3) 


Uu 


A . ~,2 i ^2 i ^2 „i,2 ri 

J\ . X -r y -r z — w — U 
B : 2x 2 + Ay 2 - w 2 = 


iJ 

/nil 1 1 Q 1 A \ 

(0|1|2|3|4) 


(0,(0, 3), 1,(1, 2), 2, (2,1), 3, (3,0), 4) 




A: x 2 + y 2 + z 2 — w 2 = 
B : 2x 2 + 4y 2 + 3z 2 - w 2 = 


[Hll]2 


iJ 

<1|2|3> 


(1,(1,2),2,(2,1),3) 


U 


A - 9t"i/ 4- _i_ ,„2 — n 

B : -x 2 + y 2 + z 2 + 2«i 2 = 


[Hll]o 


4 | 

(2) 


(2) 


\ ) 


.4 : xy + zw = 

B : -x 2 + y 2 - 2z 2 + zw+ 

2w 2 = n 


[2H]s 


5 1 

— 1 <2;;_ 2|3|2> 

(2ii + 2|3|2) 


(2,((2,1)),2,(2,1),3,(2,1),2) 
(2,((1,2)),2,(2,1),3,(2,1),2) 




A: x 2 - y 2 + z 2 + 4j/iu = 
B : -3x 2 + y 2 + z 2 = 


_6j 

(1!!_1|2|3) 


(1,((1,2)),1,(1,2),2,(2,1),3) 


• 


A : -x 2 - z 2 + 2j/«) = 
B : -3x 2 + y 2 - z 2 = 


1A 


(1,((0,3)),1,(1,2),2,(2,1),3) 


• 


4 : z 2 + z 2 + 2yw = 
B : 3x 2 + y 2 + z 2 = 


[211]! 


8 1 

(2;;_2) 


(2,((2,1)),2) 


X 


*4 : xy + zu> = 

B : 2xj/ + ?/ 2 - z 2 + w 2 = 


[22] 2 


9 1 

— \2ll_2ll_2) 
(2ll_2tt,2) 


(2,((2,1)),2,((2,1)),2) 
(2,((2,1)),2,((1,2)),2) 




A : + zu> = 

B : y 2 + 2zw + w 2 = 


[22] o 


10| 

(2) 


(2) 




A : xw + yz = 
B : xz — yw = 


[31] 2 


III 

(li«+2|3> 


(1,(((1,2))),2,(2,1),3) 


o 


A: y 2 + 2xz + w 2 =0 
B : 2yz + w 2 = 


[4]l 


12] 

(2;;;;_2) 


(2,((((2,1)))),2) 




„4 : xio + j/z = 
B : z 2 + 2?/w) = 



We sometimes also consider two simpler quadrics in the pencil spanned by A 
and B. Note that any two distinct members of the pencil have the same QSIC 
as that of A and B, and their characteristic polynomial is only different from 
that of A and B by a projective (i.e., rational linear) variable substitution. 



2.2 Open curve components 

If a connected or an irreducible component C of a QSIC is intersected by 
every plane in PR 3 , then C is called an open component; otherwise, C is called 
a closed component. A closed component curve C is compact in some affine 
realization of PR 3 ; such an affine realization is obtained by designating the 
plane at infinity to be a plane in PR 3 that does not intersect C. For example, 
a real non-degenerate conic C is closed in PR 3 . In contrast, an open component 
curve is unbounded in any affine realization of PR 3 ; a real line, for example, is 
an open curve in PR 3 . Another familiar example in PR 2 , 2D projective plane, 
is given by a nonsingular cubic curve with two connected components; it is 
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Table 2 

Classification of planar QSIC in PR 3 - Part I 



[Segre] r 
r = the # 
of real roots 


Index 
Sequence 


Signature Sequence 


Illus- 
tration 


Representative 
Quadric Pair 


[(11)"]3 


13] 

(2||2|1|2) 


(2,((1,1)),2,(1,2),1,(1,2),2) 




a 


) 


^4 : x 2 - j/ 2 + z 2 - to 2 = 
B: x 2 - 2y 2 = 


14j 

(1||3|2]3) 


(1,((1,1)),3,(2,1),2,(2,1),3) 


( 




3 


A ■ -x 2 + y 2 + z 2 + w 2 = 
B : -x 2 + 2y 2 = 


15] 

(1||1|2|3) 


(1,((0,2)),1,(1,2),2,(2,1),3) 


V X / 


^ : x 2 + y 2 + z 2 - w 2 = 
8 : x 2 + 2y 2 = 


16] 

(0||2|3|4> 
(1||3|4|3> 


(0,((0,2)), 2,(2,1), 3, (3,0),4) 
(1,((1,1)),3,(3,0),4,(3,0),3) 




^4 : x 2 + y 2 - z 2 - w 2 = 
B : x 2 + 2j/ 2 = 


[(ll)ll]l 


17| 

(1||3) 


(1,((1,1)),3) 


( 






.4 : x 2 + y 2 + 2zw = 
B : -z 2 +w 2 + 2zw = 


18] 

(2||2> 


(2,((1,1)),2) 




A : x 2 — y 2 — 2zw = 
B : -z 2 +w 2 + 2zw = 


[(lll)l]a 


19] 

(1|||2|3> 


(1,(((0,1))),2,(2,1),3) 




c 


) 


A : y 2 + z 2 -w 2 = 
B : x 2 = 


20] 

1 n Mini a \ 

(0|||3|4) 


(0,(((0,1))),3,(3,0),4) 


1 » 
1 1 

V I 

v / 


/l ■ _i_ -,2 i , ,2 — ri 

B : x 2 = 


[(21)1]2 


2l] 

<1H_|2|3) 


(1,(((1,1))),2,(2,1),3) 


( 




D 


^4 : j/ 2 - z 2 + 2zw = 
B : -x 2 + z 2 = 


22] 

(lii+|2|3> 


(1,(((0,2))),2,(2,1),3) 


[ 1 ) 


.A : j/ 2 — z 2 + 2zto = 
B : x 2 + z 2 = 



well known that one of the two components is open (i.e., intersected by every 
line in PR 2 ) and the other one is closed (i.e., not intersected by some line in 
PR 2 ). We will see that some higher order open curve components occur in 
several QSIC morphologies. 

We stress that whether a curve component is open or closed is a projective 
property, i.e., this property is not changed by a projective transformation to 
the curve. Therefore we need to consider it for classification of QSICs in PR 3 . 
In fact, the name "open component" is used here due to the lack of a more 
appropriate name, because any irreducible or connected component of a QISC 
is always "closed" in the sense that it is homomorphic to a circle. However, in 
this paper we consider the equivalence of two curve components PR 3 from the 
point of view of isotopy, i.e., homotopy of homomorphisms, as used in for knot 



theory (13 ll ). In this sense, an open curve component (i.e., intersected by every 



9 



Table 3 

Classification of planar QSIC in PR 3 - Part II 



[Segre] r 

r = the # 
of real roots 


Index 
Sequence 


Signature Sequence 


Illus- 
tration 


Representative 
Quadric Pair 


[2(H)]2 


23j 

(2!!_2||2) 


(2,((2,1)),2,((1,1)),2) 


\ / 

6 


A: 2xy - y 2 = 
B : y 2 + z 2 -w 2 = 


24J 

(1H-1II3) 


(1,((1,2)),1,((1,1)),3) 




A: 2xy - y 2 = 

B : j/ 2 - z 2 - w 2 = 


25J 

(1H+1II3) 


(1,((0,3)),1,((1,1)),3) 




A : 2xy - y 2 = 

B : j/ 2 + z 2 + w 2 = 


[(31)] 1 


26j 

(2;«_|2> 


(2,((((1,1)))),2) 




A: y 2 + 2xz - w 2 = 
8 : yz = 


27J 

(liii+|3> 


(1,((((1,1)))),3) 




.4 : j/ 2 + 2xz + w 2 = 
B : = 


[(H)(H)]2 


28 1 

(2||2||2) 


(2,((1,1)),2,((1,1)),2) 




.4 : x 2 - y 2 = 
B : z 2 - w 2 = 


291 

(0||2||4) 


(0,((0,2)),2,((2,0)),4) 




.A : x 2 + y 2 = 
B : z 2 + ™ 2 = 


30 1 

(1||1||3> 


(1,((0,2)),1,((1,1)),3) 


>c >c 


A: x 2 +y 2 = 
B: z 2 -w 2 = 


[(H)(ll)]o 


31 1 

(2) 


(2) 




„4 : xy + zw = 

B : -x 2 +y 2 - z 2 + w 2 = 


f(211)li 

L\ / J -L 


QO 1 

OZ | 

(2H-II2) 


(2,((((1,0)))),2) 




A: x 2 - y 2 + 2zw = 
B : z 2 = 


33J 

(1!!_||3> 


(1,((((1,0)))),3) 




A: x 2 + y 2 + 2zuj = 
B : z 2 = 


[(22)] i 


34j 

(2«_«_2) 


(2,((((2,0)))),2) 




A : xy + zw = 
B : j/ 2 + ui 2 = 


35J 

(2«_« + 2> 


(2,((((1,1)))),2) 




A : - zuj = 
B : j/ 2 — ui 2 = 



plane in PP 3 ) and a closed component (i.e., not intersected by some plane in 
PR 3 ) are not equivalent, because they cannot be mapped into each other by 
an isotopy of PP 3 . 
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2.3 Simultaneous block diagonalization 



When given two arbitrary quadrics, we use a projective transformation to si- 
multaneously map the two quadrics to some simpler quadrics having the same 
QSIC morphology and the same root pattern of the characteristic equation. 
Such a projective transformation is based on the standard results on simulta- 
neous block diagonalization of two real symmetric matrices (1251 ; l49l; |40|) . which 
will be reviewed below. 

Definition 1: Let A and B be two real symmetric matrices with A being 
nonsingular. Then A and B are called a nonsingular pair of real symmetric 
(r.s.) matrices. 



Definition 2: A square matrix of the form 



M 



A e 



V 



e 
A 



kxk 



is called a Jordan block of type I if A G R and e = 1 for k > 2 or M — (A) 
with A G R for k = 1; M is called a Jordan block of type II if 




a,b e R, b and e 




for k > 4 or 



M 




for k = 2, with a,b e R, b ^ 0. 

Definition 3: Let Ji,...,J^ be all the Jordan blocks (of type I or type II) 
associated with the same eigenvalue X of a real matrix A. Then 

C = C(A) = diag(J 1 ,...,J fc ), 

where dim(Jj) > dim(Jj + i) ; is called the full chain of Jordan blocks or full 
Jordan chain of length k associated with A. 

Definition 4: If Xi,...,\k are all distinct eigenvalues of a real matrix A, with 
only one being listed for each pair of complex conjugate eigenvalues, then the 
real Jordan normal form of A is J=di&g(C(Xi),...,C(Xk)). 
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Recall that two square matrices C and D are congruent if there exists a non- 
singular matrix Q such that C = Q T DQ; we also say that C and D are related 
by a congruence transformation, which amounts to a change of projective co- 
ordinates. 

Theorem 1 (Quadric Pair Canonical Form) 

Let A and B be a nonsingular pair of real symmetric matrices of size n. Sup- 
pose that A~ X B has real Jordan normal form diag(Ji, ...J r , J r +i, ■■■J m ), where 
Ji,...J r are Jordan blocks of type I corresponding to the real eigenvalues of 
A~ X B and J r+ \, ...Jm ore Jordan blocks of type II corresponding to the com- 
plex eigenvalues of A~ X B. Then the following properties hold: 

(1) A and B are simultaneously congruent by a real congruence transforma- 
tion to 

diag(£i-Bi, ...e r E r , E r+ %, ...E m ) 

and 

diag(ei-Ei J\, ...e r E r J r , E r+ iJ r+ i, ...E m J m ), 
respectively, where £j = ±1 and the Ei are of the form 

^0.0 1^ 



v io.o, 

of the same size as J i} i — 1,2, .., m. The signs of Si are unique for each 
set of indices i that are associated with a set of identical Jordan blocks Ji 
of type I. 

(2) The characteristic polynomial of A~ X B and det(XA — B) have the same 
roots Xj with the same multiplicities ji. 

(3) The sum of the sizes of the Jordan blocks corresponding to a real root A, 
is the multiplicity 7, if \ is real or twice this multiplicity if \ is complex. 
The number of the corresponding blocks is pi = n — rank(Ajv4 — B), and 
1 < Pi < Ji- 



As detailed in the review article (16), this result has a long story It was proved 
in (jiil ) for non-degenerate pencils, and them further extended and rediscovered 



several times. See (135 



ate pencils, ana tnem : 



In order to apply Theorem 1, we need to ensure that the matrix A is nonsin- 
gular. Since we assume that /(A) = det(Av4 — B) does not vanish identically, 
XA — B is nonsingular for infinitely many values of A. Therefore, given two 
quadrics A : X 7 AX = and B : X T BX = 0, we may assume that A is 
nonsingular; for otherwise we may replace A by another nonsingular matrix 
A such that A : X T AX = and E have the same QSIC as that of A and B. 
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2.4 Index sequences 



Signature and index: Any n x n real symmetric matrix D is congruent to 
a unique diagonal form D' = diag(ij, —Ij, 0^). The signature, or inertia, of D 
is (<7+, <r_, cr ) = («, j, fc). The index of .D is defined as index(D) = i. 

Index function: The index function of a quadric pencil XA — B is defined as 

Id(A) = index(XA-B), A G PR. 

Since A and B are matrices of order 4 in our discussion, i.e., n = 4, we have 
Id(A) G {0,1,2,3,4}. Note that Id(A) has a constant value in the interval 
between any two consecutive real roots of /(A) = 0. The index function may 
have a jump across a real root of /(A) = 0, depending on the nature of 
the root. The index function is also defined for A = oo and — oo. We have 
Id(-oo) + Id(+oo) = rank(A). 

Eigenvalue Curve: We consider the real eigenvalues of the pencil \A — B, 
defined by the equation 

C(A, u) = det(\A - B - u I) = 0. 

We are going to see that the QSIC of a pencil (A, B) can be characterized 
by the geometry of the planar curve C defined by the equation C(A, u) = 0. 
This curve C is defined by a polynomial whose total and partial degree in 
either A or u is 4. Since a 4 x 4 symmetric matrix has 4 real eigenvalues, for 
any A G M, the number of real roots C(A, u) = in u is 4 (counted with 
multiplicities). Consequently, there are 4 A-monotone branches of C. For any 
fixed A G M, the number of points of C not on the A-axis, i.e., with u ^ 0, is 
the rank of the quadratic form XA — B; the number of points of C above the 
A-axis and the number of points of C below the A-axis determine the signature 
of (A A — B). Figure 1 shows the eigenvalue curve of the pencil of quadrics 
(y 2 + 2xz+ 1, 2yz + l). 

Index sequence: Let Aj , j = 1,2,..., r, be all distinct real roots of /(A) = 
in increasing order. Let k = 1, 2, . . . , r — 1, be any real numbers separating 
the Xj, i.e., 

— oo < A x < fii < X 2 < . . . /i r -i < A r < oo. 

Denote Sj = ld(/ij), j = 1, 2, . . . , r — 1. Denote s = Id(— oo) and s r = Id(oo). 
Then the index sequence of A and B is defined as 

( s T s l T • • • T s r-l T s r ), 
where j stands for a real root, single or multiple, of /(A) = 0. 
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Fig. 1. The eigenvalue curve of the pencil of the quadrics (y 2 + 2 x z + 1, 2y z + 1). 

To distinguish different types of multiplicity of a real root, we use | to denote 
a real root associated with a 1 x 1 Jordan block, and use I for p consecutive 
times to denote a real root associated with apxp Jordan block. For example, 
a real root with Segre characteristic [11] will be denoted by || in place of an j 
in the index sequence, and a real root with the Segre characteristic [21] will be 
denoted by ll\ in place of an |. When the Segre characteristic is (22), we use 
Ull to distinguish it from Ull, which has the Segre characteristic [4]. Supposing 
that Ao is a real zero of /(A) with a Jordan block of size kxk,we use I • ■ • l + or 
I - ■ - l- to indicate that the corresponding sign Si of the block in the Quadric 
Pair Canonical Form is + or — . 

Since A is a projective parameter, a projective transformation A' = (aA + 
b)/(c\ + d) does not change the pencil but may change the index sequence 
of the pencil. On the other hand, thinking of the projective real line of A 
as a circle topologically, such a transformation induces either a rotation or 
a reversal of order of the index sequence of the pencil. Therefore we need 
to define an equivalence relation of all index sequences of a quadric pencil 
under projective transformations of A. In addition, replacing A and B by —A 
and —B changes each index Sj to rank(AA — B) — s$ but essentially does not 
change the pencil XA — B. Note that the above replacement changes of the 
sign associated with a Jordan block of a root; for instance, if the quadrics A 
and B have the index sequence (2?L2|3|2), then —A and —B have the index 
sequence (2» + 2|l|2). 

We choose a representative in an equivalence class such that A is nonsingular; 
therefore, oo is not a root of /(A) = and s + s r = 4. Taking these observa- 
tions and conventions into consideration and denoting the equivalence relation 
by ~, this equivalence of index sequences is then defined by the following three 
rules: 

1) Rotation equivalence: 
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(s Q t si t • • • T s r _i t s r ) ~ (4 - s r _i t s t -si t • • • T s r _i), (2) 

(S T Si T • • • T Sr-l T Sr) ~ («1 T «2 T • • • T S r t 4 - Si). 

2) Reversal equivalence: 

(S Q t Si t • • • T Sr-l t S r ) ~ (S r t S r _i t • • • T Si t S ). (3) 

3) Complement equivalence: 

(S T Sl T • • • T Sr-l T S r ) ~ (4 - S T 4 - Si t ... T 4 - S r -1 | 4 - S r ). (4) 

£.5 Signature variation 



In this section we analyze the behavior of the eigenvalues of the pencil H(X) = 
XA — B, near the roots of /(A) = det(H(X)) = 0. This analysis amounts 
to analyzing the eigenvalue curves at a real root of /(A), and is needed for 
computing the jump of the index function at the real root. 

Consider a transformation H'(\) = P T H(X)P of H(X), where P is an invert- 
ible matrix. First, we compare the behavior of the eigenvalues of H'(\) and 
H(X). For any real symmetric matrix Q of size n, we denote by Pk{Q) the 
k th real eigenvalue of Q, so that pi(Q) ^P2(Q) ^ • ■ • ^ Pn{Q)- Using the 
Cour ant- Fischer Maximin Theorem (see (1141 ) p. 403), we have the following 
result: 

p fc (QVi(P) 2 < Pk (P T QP) ^ p k {Q)a n {Pf, (5) 
where <Ji(P) (resp. a n (P)) is the smallest (resp. largest) singular value of P. 

Proposition 1 Let P be an invertible matrix and H'(X) = P T H(X)P. If 
p k (H(X)) = aX^(l + o(A)) with a ^ 0, then p k (H'(X)) = a'X^(l + o(A)) with 
sign(a) = sign(a'). 



Proof As the eigenvalue pk(H'(X)) has a Puiseux expansion yj; |42j) near 
A = of the form p k (H'(X)) = p' + a'A^'(l + o(A)) with p', s'el and // G Q, 
we deduce from the inequalities (5) that p' — 0, p! — p and sign(a') = sign(a). 

Proposition 1 allows us to deduce the behavior of the eigenvalues of the pencil 
H(X), from its normal form. Indeed, by Theorem 1, H(X) is equivalent to 

D(X) = di ag (e 1 E 1 (XI 1 -J 1 (X 1 )),e 2 E 2 (XI 2 -J 2 (X 2 )), . . . ,e r E r (XI r -J r (X r )),D'(X)), 

i (6) 

where Ii is the identity matrix of the same size as that of the Jordan block 
Jj(Aj) of eigenvalue A i; and det(D'(A)) has no real roots. Let us denote by 
N k (X,p,e) = eE k (XI k — J k (p)) a block of the preceding form, where k is the 
size of the corresponding matrices. Then we have the following property: 
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Proposition 2 The eigenvalue branch p(X) corresponding to AT fc (A, po,e) which 
vanishes at X = p is of the form 

p = ev k {l + o{v)) 

where X = p + v . 

Proof By an explicit expansion of the determinant jV(A, u) = det(iVfc(A, p , e) — 
ulk) and denoting v = X — po, we obtain 

ff(\,u)=ff(v,u) = {-l) k u k + --- + {-e) k -\-l) {J ^^ +1 u+e\-l)^u\ 

The vertices of the lower envelop of the Newton polygon of Af(v, u) in the 
(u, v) -monomial space are the points (k, 0), (1, 0), (0, k). By Newton's theorem 
(see (jH) p. 89), the Puiseux expansion of the root branch which vanishes near 
Po is of the form 

p = ei/ fc (l + o(i/)), 

which completes the proof. 

According to Proposition 1, if the pencil H(X) is equivalent to (6), then near 
each root Aj, the eigenvalue branches approaching are of the form Ei(X — 
Aj) fci (l+o(A — Aj)), where ki is the size of a block of the Quadric Pair Canonical 
Form (6) of the eigenvalue Aj and is the corresponding sign. 

Index Jump: The preceding analysis explains how the index function can 
change around the real roots of /(A) = 0. Let a be a real root of /(A) = 0. 
Let a_ and a + be values sufficiently close to a, with a_ < a and a + > a. 
Then the index jumps of Id(A) at a are denoted as 

A" (a) = Id(a) - Id(a_), A+(a) = Id(a+) - Id(a), 
A(a) = Id(a+) - Id(a_) = A" (a) + A+(a). 

We denote by Af (a) the changes of signature functions of the blocks iV^A, Aj, Si) 
at a. Clearly we have 

k k 

A ± (a) = EAf(a), A(«) = ^A(«). (7) 

i=l i=l 

Let us describe each Af (a) separately. For any a £ R, we denote a + = 
max(a, 0) and a~ = min(a, 0). Note that a + + a~ = a. 

(1) Jordan block of size 1 x 1: In this case, clearly, we have the following 

signature sequence (e~ , (0,0), ef) and the jumps are A~(A«) = e~, A + (Aj) = 
ef and A(Aj) = £j. 
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(2) Jordan block of size 2x2: 



N k .(X, Ai,ei) 



A- A; 
A - A,- -1 



In this case the corresponding eigenvalue branch vanishing at Aj is equivalent 
to £j(A — Aj) 2 ; therefore its sign is the same before and after Aj. There is one 
positive eigenvalue and one negative eigenvalue before and after Aj. If £j > 0, 
we have a positive eigenvalue branch which goes to at Aj; otherwise, we have 
a negative one. Thus, the signature sequence of N k .(\, Aj, £j) is (1, (1 — ef, 1 + 



£j ), 1) and the jumps are A (Aj) 



-e+ A+(A i )=ef andA(A i ) = 0. 



(3) Jordan block of size 3x3: Since 



N ki (\, A,, Si) — Si 



A - \ 
A-Ai -1 
A - A,- -1 



The corresponding eigenvalue branch is equivalent to £j(A — Aj) 3 , whose sign 
changes before and after Aj. If £j > 0, the signature of N k .(X, Aj, £j) is (1,2) 
before Aj and (2, 1) after. If £j < 0, we exchange the order of the two signatures. 
Thus, we have the signature sequence (ef — 2e~ , (l,l),2ef — e^) — (1 — 
£-, (1, 1), 1 + ef) and A~(Aj) = e~ , A+(A,) = ef and A(A,) = e l . 

(4) Jordan block of size 4x4: Using a similar argument, we can show that 
there are two positive eigenvalues and two negative eigenvalues before and 
after Aj and the eigenvalue curve approaching zero has the form £j(A — Aj) 4 . 
Thus, the signature sequence of N k .(\,\ i: Ei) is (2,(2 — ef,2 + el), 2) and 
A-(Aj) = -£+,A+(Aj) = £+,A(Aj) = 0. 

To summarize, taking into account the sign £j = ±1, we have Aj(a) = £j if Jj 
has the size 1 x 1 or 3 x 3, and Aj(a) = if Jj has the size 2 x 2 or 4 x 4. The 
rank of H(\) drops by 1 at A = Aj for each block of the form N k .(\\ h Ei). 
Thus, the signature of iJ(Aj) can be deduced directly from its index Id(Aj) 
and the number of Jordan blocks with eigenvalue Aj. 

The above rules can be used to decide the permissible index jumps of Id(A) 
at a real root of /(A) = 0, through Eqn. (7) and the signature of H(\i). In 
particular, in the case of a simple root Aj of /(A) = 0, the sign £j in the Quadric 
Pair Canonical Form can be deduced directly from the index before and after 
the root. For instance, an index sequence of the form (1|2|1|2|3) corresponds 
to a sequence of signs e\ = +1, £2 = — 1, £3 = +1, £4 = +1, and the signatures 
at the roots are (1,2), (1, 2), (1, 2), (2, 1), respectively. 
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Signature sequence: The previous analysis allows us to completely deter- 
mine the signature sequence of the pencil H(X) = XA — B, from its Quadric 
Pair Canonical Form. For most of the cases, this signature sequence is, as we 
will see, a characterization of the QSIC. A signature sequence is defined as 

(s Q , (• ■ ■ (pi.ni) • • • ), sx, ■ ■ ■ s r _i, (• ■ ■ (p r , n r ) ■■■), s r ), 

where Sj is the index of H(X) between two consecutive real roots of /(A) = 0, 
(pi,rii) is the signature of H(\i) at a root Aj and the number of parentheses 
is the multiplicity of Aj. Note that Pi + = rank(AjA — B). 

The advantage of using the signature sequence over using the index sequence 
is that we just need to compute the multiplicity of a real root and determine 
the signature of A A — B at the root; this is a far simpler computation than 
computing the Jordan block size, which is the information required by the 
index sequence. Conversion from an index sequence to the corresponding sig- 
nature sequence is straightforward. For a given pair of quadrics, the signature 
sequence can be computed easily using only rational arithmetic as described in 
Section 2.6. Similar equivalence rules to those for index sequences apply to sig- 
nature sequences as well. The signature sequences of all 35 QSIC morphologies 
are listed in the third column of Tables 1, 2 and 3. 

2.6 Effective issues 

Now we discuss how to use rational arithmetic to compute the signature se- 
quence for classifying the QISC morphology of a given pair of quadrics. Con- 
sider the polynomial 

C(A, u) = det{XA — B — ul) = w 4 + c 3 (A)w 3 + c 2 (A)w 2 + a(X)u + c (A). 

The values where the signature changes are defined by C(A, 0) = c (A) = 
/(A) = 0. For a fixed A, the rank of the corresponding quadratic form is 
the number of non-zero roots of C(A, u) = 0. For any fixed A, the number 
of real roots in u, counted with multiplicity, is 4. The signature of XA — B 
is determined by the rank of XA — B and the number of positive roots of 
C(A, u) = in u. In the case where the number of real roots equals the degree 
of the polynomial, the Descartes rule gives an exact counting of the number 
of positive roots p), and we have the following property: 

Theorem 2 For any X E R, 

• the number of positive eigenvalues of XA — B is the number of sign variations 
o/[l,C3(A),c 2 (A), Cl (A),co(A)]. 

• the number of negative eigenvalues of XA—B is the number of sign variations 
o/[l,- C 3(A),c 2 (A),- Cl (A),co(A)]. 



18 



Computing the signature XA — B for A G Q is straightforward. Computing its 
signature at a root of C(A,0) = /(A) = can also be performed using only- 
rational arithmetic. According to the previous propositions, this reduces to 
evaluating the sign of Cj(A), z = 1 ... 3. This problem can be transformed into 
rational computation as follows. First, we represent a root a of /(A) = by 

• the square-free part p(X) of /(A) = and 

• an isolating interval [a, b] with a, b G Q such that a is the only root of p(\) 
in [a, 6]. 

Isolating intervals can be obtained efficiently in several ways (see, for instance, 
(231)). They can even be pre-computed in the case of polynomials of degree 4 
(llOl ). In order to compute the sign of a polynomial g at a root a of /(A) = 
0, we use subresultant (or Sturm-Habicht) sequences. We recall briefly the 
construction here and refer to (|2j) for more details. 

Given two polynomials /(A) and g(X) G A[A], where A is the ring of coeffi- 
cients, we compute the sub-resultant sequence in A, defined in terms of the 
minors of the Sylvester resultant matrix of /(A) and f'(X)g(X). This yields a 
sequence of polynomials R(A) = [R (X), Ri(X), . . . , Rn(X)] with Ri(X) G A [A], 
whose coefficients are in the same ring A. 

In our case, we take A = Z. For any a G R, we denote by V/ i9 (a) the number 
of sign variation of R(a). Then we have the following property (|2j): 

Theorem 3 

V f,g( a ) ~ V f,g( b ) = #{« £ [°> b] root of /(A) = where g(a) > 0} - 
#{a G [a, b] root of /(A) = where g(a) < 0}. 



In particular, if the interval [a, b] is an isolating interval for a root a of Cq(X) = 
0, then Vf t g(a) — Vf t g(b) gives the sign of g(a). Taking g(X) to be the coefficients 
q(A) in Theorem 2, this method allows us to exactly compute the signature 
of a A — B, using only rational arithmetic. 

Efficient implementations of the algorithms presented here are available in the 
library synap£] and have been applied to classifying QSIC morphologies, 
based on the signature sequences derived in this paper. 



http : //www-sop . inria.fr/galaad/software/ synaps/ 
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2. 7 List of QSIC morphologies 

All 35 different morphologies of QSIC are listed in Tables 1 through 3. In the 
first column are the Segre characteristics with the subscript indicating the 
number of real roots, not counting multiplicities. The index sequences and 
signature sequences are given in the second column and the third column, 
respectively. Here, only one representative is given for each equivalence class 
associated with the corresponding QSIC morphology; in several cases, there 
are two equivalence classes associated with one QSIC morphology. The nu- 
meral label for each case, from 1 to 35, is given at the left upper corner of 
each entry in the second column. These labels are referred to in subsequent 
theorems establishing the relation between the index sequence and the QSIC 
morphology. Cases 4, 10 and 31 share the same index sequence (2), thus also 
the same signature sequence (2). Additional simple conditions based on mini- 
mal polynomials for distinguishing these three cases are presented in Section 6. 
Two different index sequences in cases 26 and 34 correspond to the same sig- 
nature sequences; the discrimination of these two cases is also discussed in 
Section 6. In the illustration of each QSIC morphology in column four, a solid 
line or curve stands for a real component and a dashed one depicts an imag- 
inary component. A solid dot indicates a real singular point, which in many 
cases is a real intersection point of two or more components of a QSIC. An 
open or closed component is drawn as such in the illustration. Note that, 
in addition to topological properties, we also take algebraic properties into 
consideration in defining morphologically different types. For example, a non- 
singular QSIC may be vacuous in PR 3 , so is a QSIC consisting two imaginary 
conies; these two QSICs are defined to be morphologically different since the 
former is irreducible algebraically but the latter is not. 



3 Classifications of nonsingular QSIC 

3.1 [111 1)4. • /(A) = has four distinct real roots 

Theorem 4 Given two quadrics A: X T AX = and B: X T BX = 0, if their 
characteristic equation /(A) = has four distinct real roots, then the only 
possible index sequences are (1|2|1|2|3) and (0|1|2|3|4). Furthermore, 

(1) (Case 1, Table 1) when the index sequence is (1|2|1|2|3) 7 the QSIC has 
two closed components; 

(2) (Case 2, Table 1) when the index sequence is (0|1|2|3|4) ; the QSIC is 
vacuous in PR 3 . 
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Proof Let Aj, i = 1,2,3,4, be the four distinct real roots of /(A) = 0. By 
Theorem 1, A and B are simultaneously congruent to 



A = diag(ei, e 2 , e 3 , e 4 ), and B = diag(ejAi, e 2 A2, £3X3, £4A 4 ) 



where Si = ±1, i = 1,2,3,4. Without loss of generality, we suppose that 
Ai < A2 < A3 < A4; this permutation of the diagonal elements can be achieved 
by a further congruence transformation to A and B. 

Clearly, the only possible index sequences are (up to the equivalence rules of 
Section 2.4) (1|2|1|2|3) and (0|1|2|3|4). Since a pencil with the second index 
sequence (0|1|2|3|4) contains a positive definite or negative definite quadric, 
i.e., with the index being 4 or 0, we deduce that the intersection curve is empty 
in that case. 

For the first index sequence (1|2|1|2|3), according to Section 2.5, the sign 
sequence in the corresponding Quadric Pair Canonical Form is (ei = I, £2 = 
— 1, £3 = 1, £4 = 1). Setting A to A' and B — X4A to B', we obtain 

A' = diag(l,-l,l,l), 



B' = diag ((Ai - A 4 ), -(A 2 - A 4 ), (A 3 - A 4 ), ) . 

Consider the affine realization of PR 3 by making y = the plane at infinity. 
Then A' is a sphere, which intersects the x-z plane in a unit circle, while the 
quadric B' is an elliptic cylinder with the w-axis being its central direction, 
which intersects the x-z plane in an ellipse, since \ < A 4 , i = 1, 2, 3. Clearly, if 
one of the ellipse's semi-axes is smaller than 1 or both are smaller than 1, the 
QSIC of A' and B' has two oval branches (see the left and middle configurations 
in Figure 2). If both of the ellipse's semi-axes are greater than 1, A' and B' 
have no real intersection points (see the right configuration in Figure 2). We 



recall the following result from (jilt l39f ) : Two quadrics A : X 7 AX = and 
B : X T BX = in PM 3 has no real points if and only if XqA — B is positive 
definite or negative definite for some real number X . It implies that the index 
sequence of the pencil cannot be (1|2|1|2|3). This is a contradiction. Hence, 
the QSIC has two ovals. 

Note that none of the semi- axes can be of length 1, since /(A) = is assumed to 
have no multiple roots. We deduce that the QSIC has two closed components 
when the index sequence is (1|2|1|2|3) and is empty when the index sequence 
is (0|1|2|3|4). This completes the proof of Theorem 4. 
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Fig. 2. Three cases of an elliptic cylinder intersecting with a unit sphere and their 
corresponding cross sections in the x-z plane. 

3.2 [1111] 2 .' /(A) = has two distinct real roots and a pair of complex 
conjugate roots 




Theorem 5 (Case 3, Table 1) If f(X) = has two distinct real roots and one 
pair of complex conjugate roots, then the index sequence of the pencil XA — B 
is (1|2|3) ; and the QSIC comprises exactly one closed component in PR 3 . 

Proof Wlog, we assume A is nonsingular. Suppose that /(A) = has two real 
roots Ai 7^ A2 and two complex conjugate roots A34 = a±bi. First, it is easy 
to see that the only index sequence possible is (1|2|3). We may suppose that 
A3 ; 4 = ±i; this can be done by setting (B — aA)/b to B. By Theorem 1, A and 
B are congruent to 



A' 



diag(£d,£i,£ 2 ) = dia S II I : 'i 

1 



B ' = (Kj) = diag(E 1 J 1 ,e 1 \ 1 ,e 2 \ 2 ) 
As the index sequence is (1|2|3), we have S\ 
two cases: (1) AiA 2 7^ and (2) AiA 2 = 0. 




-1, 1, £]Ai, £ 2 A 2 ). 
1. Next we need consider 



Case 1 (AiA 2 7^ 0): By a variable transformation A' = —A if necessary, we 
may assume that at least one of Ai and A 2 is positive. Then we denote Ai > 
and A 2 < if only one of them is positive or denote A 2 > \\ > if both are 
positive. It follows that ^ < 1. We then set XiA' — B' to A' and use a further 
simultaneous congruence transformation to scale the diagonal elements of B' 
into ±1. For simplicity of notation, we use the same symbols A' and B' for 
the resulting matrices and obtain 



A' 



131 



1 Xi 
Ax -1 



V 



where /5 2 = A 2 /|A 2 | = ±1. 



\ 



B' = {b' t 



hi 
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If /?2 = 1, we swap 644 and b' lx , as well as a' A A and a' l x , to obtain 



A' 





\ 




(l 


\ 




-1 X, 


, B' = 


1 











1 




V 


Ai 1 ) 




V 


- 1 ) 



Or, if fa — — 1, we swap 6 4j4 and & 22 , as well as a' A i and a 22 , to obtain 



A' = 



Ai 





(-1 


\ 






-1 


, B' = 










1 


J 


\ 





Note that permuting diagonal elements can be achieved by a congruence trans- 
formation. Hence, whether /3 2 = 1 or /3 2 = — 1, after a proper simultaneous 
congruence transformation, i3' is the unit sphere or a one-sheet hyperboloid 
with the z-axis as its central axis. Since ^ < 1, and a 2 2 have the same 
sign. Therefore, A! is an elliptic cylinder parallel to the z-axis. Due to the 
symmetry of B' and A 1 about the x-y plane, we just need to analyze the rela- 
tionship between the two conic sections in which A' and £>' intersect with the 
x-y plane. 

The quadric £>' intersects the x-y plane in the unit circle x 2 + y 2 = w 2 , and 
A' intersects the x-y plane in the ellipse 



a (y 

a 2 b 2 



cw) 



= w 



when — 1, or in the ellipse 

(x + cw) 2 y 2 



b 2 



a 1 



when /3 2 — — 1- Here a — -1 / t|~^j ' ^ = \/l + A?, and c = Ai. 



In both cases of @ 2 — =tl, the center of the ellipse shifts from the origin 
(along the x direction or y direction) by the distance |Ai|, and the length 

of the ellipse's semi-axis in the shift direction is b = \Jl + A 2 . Then it is 
straightforward to verify that one of the ellipse's extreme points of this axis 
is inside the unit circle, while the other is outside the unit circle. (See Figure 
3 for the case of ^2 = —1.) In this case the QSIC of A' and B' has one closed 
component in PM 3 (see Figure 4). 
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Fig. 3. The cross-sections of an elliptic Fig. 4. The intersection curve referred to 
cylinder and a hyperboloid with one sheet in Figure 3 
in the x-y plane. 

Case 2 (A1A2 = 0): Wlog, we may suppose that Ai = and A2 7^ 0. Then, by 
Theorem 1, noting that e\ — £2 = 1, A and B are congruent to 



A' 



foi 

1 

V 



\ 



B' 



I 



-1 
1 



\ 



V 



A 



2/ 



First set A' — (1/X 2 )B' to be A'. Then we use a congruence transformation 
to make the diagonal elements of B' become ±1 and apply the same trans- 
formation to A'. Denoting the resulting matrices again using A' and B', we 
obtain 



A' 



Ik ' 
1 -A 



\ 



( 



-1 
1 



We swap b' 4 A and U X1 , as well as a' A 4 and a' tl , by a simultaneous congruence 
transformation to obtain 



A' 







\ 



j_ 

A 2 



1 i 



B 1 



A 2 / 



1 



Thus, B' is a cylinder with the z-axis as its central axis, and A! is either an 
elliptic cylinder or a hyperbolic cylinder, depending on the sign of A2, and A 1 
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is parallel to the y-axis. The equation of A' is 

(y - cwf z 2 2 

2 ± Ta = w ' 

a 2 ¥ 

where a = \J\ + b = -j^f-, c = X 2 . The cylinder A' shifts from the origin 
by the distance | A2 1 along the x-axis or the y-axis, and the length of its semi- 
axis in the shift direction is J 1 + A| . Clearly, in this case, the QSIC of the 
cylinders A' and £>' has exactly one closed component in PR 3 . (See Figure 5.) 
This completes the proof. 




Fig. 5. The intersection of a circular cylinder with a hyperbolic cylinder or an elliptic 
cylinder. 



3.3 [ 1 1 1 1 ] .' /(A) = has two distinct pairs of complex conjugate roots 



Theorem 6 (Case 4, Table 1) If f(A) = has two distinct pairs of complex 
conjugate roots, then the Segre characteristic is [1111] and the index sequence 
is (2). In this case the QSIC comprises two open components in PR 3 . 

Proof Suppose that /(A) = has the roots a ± bi and c ± di. First, it is 
easy to see that the index sequence is (2). By setting (B — cA)/d to be B, we 
transform conjugate roots c ± di to ±z. Therefore, we suppose that /(A) = 
has the roots a±bi and ±z. Furthermore, we may suppose that A and B form 
a nonsingular pair of real symmetric matrices. Then, by Theorem 1, A and B 
have the following canonical forms 



A' = diag 




and B' = diag 



V 



—6 a 
a b 



Here, a 7^ or b 7^ ±1, since the roots a±bi are distinct from ±z. Also, b 7^ 
since a±bi are imaginary. Wlog, we may assume b > 0. In the following we will 
derive a parameterization of the QSIC from which the topological information 
about the QSIC can be deduced. The quadric A' : X T A'X = is a hyperbolic 
paraboloid and can therefore be parameterized by r(u,v) = g(w) + h(u)v 
where 

g (u) = (- M ,0,0,1) T and h(u) = (0, 1, u, 0) T . 
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Substituting r(u, v) into X B'X = yields 



(u) T B'h(u) ± ^s(u 
h{u) T B'h{u) 



where 



s(u) = [ g (u) T B'h(u)] 2 - [( g (u) T B' g (u))(h(u) T B'h(u))] 
= -bu 4 + (a 2 + b 2 + l)u 2 - b. 

Substituting (8) into r(w, v ) yields the following parameterization of the QSIC, 



p(u) 



bu 3 — u, — ( au ± J s(u)j , — u ( au ± \J s{u) j , 1 — feu 



(9) 



Since p(w) is a real point only when s(u) > 0, we are going to identify the 
intervals in which s(u) > holds. We will first show that s(u) = always has 
four distinct real roots. The equation 

s( u ) = -bu 4 + (a 2 + b 2 + l)u 2 -6 = 

is a quadratic equation in u 2 with discriminant 

A = (a 2 + b 2 + l) 2 - 46 2 = a 2 (a 2 + 26 2 + 2) + (b 2 - l) 2 > 0, 

since o ^ or 6 ^ ±1. Therefore the two real solutions of u 2 are 

2 (a 2 + 6 2 + 1) ± v/A 



2fe 



(10) 



Since A 



> 2 + fe 2 + l) 2 - 46 2 and b ^ 0, we have (a 2 + 6 2 + 1) > v^. It 
follows that the numerator and denominator in (10) are positive; recall that 
b > is assumed. Then we get the four real solutions for s(u) = from (10), 
denoted by ±u + and ±w_, with u + > w_ > 0. 





Fig. 6. The graph of s(u). 



Fig. 7. The case of the QSIC has two 
affinely infinite components. 



Define two intervals I\ = [u-,u + ], I 2 = [—u + ,—uJ\. Since s(0) = — b < 0, 
we have s{u) > for u E I\ an d s{u) < for the other values of u. (See 



26 



Figure 6 for the graph of s(u).) This implies that the QSIC, given by p(u), 
has two connected components, denoted by V\ and V 2 , corresponding to the 
intervals I\ and l 2 . V\ is defined by p(u) over the interval I±, and V 2 is defined 
by p(u) over the interval I 2 . 

Next we are going to show that the two components Vi and V 2 are open curves 
in PR 3 . Since V\ and V2 have the same parametric expression p(u) but over 
different intervals, we will only analyze the component V\, the analysis for 
V2 is similar. The key idea of the proof is to show that Vi has exactly one 
intersection point with the plane w — 0; it will then follows that V\ has an 
intersection with every plane in 



Consider the affine realization AIR 3 of PR 3 by making the plane w = the 
plane at infinity. The w-coordinate component of p(u) is w(u) = 1 — bu 2 , 
which has two zeros ui = 1/y/b and u 2 = —1/y/b, and it is straightforward 
to verify that U\ = 1/y/b G I\ and u 2 = —1/y/b £ I 2 . Therefore we will only 
consider the two points p(ui) (i.e., with ± in Eqn. (9)) on the component V\. 
Let qo(w) and qi(u) denote the two "branches" of p(u) corresponding to ± in 



front of Js(u) in Eqn. (9). Then 



q (w) = ^6-u 3 — u, — (au + \J s(u)^J , — u (au + \J s(u)^j , 1 — bu 2 ^j 

qi(-u) = (bu 3 — u, — (au — \J s(u)^j , — u (au — \J s(u)^j , 1 — bu 2 ^j . (11) 



There are now three cases to consider: (i) a = 0; (ii) a > 0; and (in) a < 0. 
First consider the case (i) a = 0. In this case, 

s( u ) = -bu 4 + (b 2 + l)u 2 - b = (u 2 - b)(l - bu 2 ). 



It follows from Eqn. (11), after dropping a common factor \J\1 — bu 2 \, that 



^\l-bu 2 \, -^\u 2 -b\, - U yJ\u 2 -b\, v/|1-M) T , 



qoW = -u 



q_i(u) = —uJ\l — bu 2 \, J\u 2 — b\,uJ\u 2 — b\,J\l — bu 2 



Note that the two ends of I\ are 1 / y/b and y/b when a = 0. It is easy to verify 
that, when u = y/b, the two branches q,o(u) and qi(u) are joined together at 
the finite point 

q (Vb) = qi (Vb) = (-^/b\l-b 2 \, 0, 0, ^\l-b 2 \) T . 
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Next consider the behavior of qo(tt) and qi(u) when u — > Ui — 1/Vb, which 
is the other end of I\. Since 



q M = (0, - 6|, - v^- 1 - 6|/6, 0) J 

and 

qaOix) = (0, Vl^-fel, y/\b~i -b\/b, 0) T , 
qo(«i) and qi(tii) represent the same point at infinity. 

Denote qj(w) = (xi,Ui, Zi,Wi) T , i = 0, 1. To study the asymptotic behavior of 
the QSIC, let us consider the limit of the affine coordinates of qo(u) and qi(u), 
i.e., (x /wo,yo/w ,zo/w ) T and (x 1 /wi,y 1 /w 1 , z 1 /w 1 ) T , as u -> ui = 
Clearly 

lim = l im £iM = _J_ 

u— ui W q(u) u ^ u i W 1 (u) y/J} 7 

,. 2/o(«) r 

lim — — = — lim — — = — oo, 

and 

r *o(") i- z i( u ) 

lim — — — = — lim — — — = — oo. 

u ~> u i w (u) Wi(u) 
Therefore, the component curve V\ comprises one connected component and 
its two ends extend to infinity in opposite directions in AM 3 , with its asymptote 
line being the intersection line of the two planes x+U\W = and — U\y + z = 0. 
Clearly, Vi is intersected by every plane in PR 3 . Hence, V\ is open in PR 3 . 

Now we consider case (ii): a > 0. In this case, ui = l/Vb £ l\ = 
Clearly, the two parts of Vi defined by qo(w) and qi(u) are joined at the two 
finite points qo(u-) — <?i( M -) an d <?o( M +) — <?i( M +) We will show that, when 
u — Ui — 1/Vb, qo(ui) gives the on/?/ infinite point on Vi. Since, from Eqn. 
(11), 

q (ui) = (0,-2a/A -2a/6,0) T , 
qo(ui) is a point at infinity. To study how qo(w) approach the infinite point 
qo(ui), let us consider the lim u _^ Ul qo(u) and lim u ^ Ul+ qo(w) when u ap- 
proaches Mi from different sides. Using the affine coordinates of qo(w), we 
have 

lim ^1 = - lim = 
lim 44 = - lim 44 = +-, 



and 



lim — — = — lim — — = +oo. 

u^m- Wo(U) u^m + Wq{U) 



Thus, the component Vi extends to infinity in opposite directions, with its 
asymptote being the intersection line of the two planes x + u±w = and 
—uiy + -2 = 0. 
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Note that all other points of Vi are obviously finite, except for the point qi(iti) 
whose w component is zero. But we will show that qi(wi) is, in fact, also a 
finite point. Denote g(u) = \Js(u). Expanding g{u) at u — U\ by the Taylor 
formula yields 

g(u) = + g'iu^iu - m) + o(u - ui), 

where o(u — U\) is a term whose order is higher than u — U\ when u — > u\. 
Plugging the above g(u) in qi(tt) yields 



qi(u) = (bu 3 — u, — (au — g(u)), —u(au — g(u)), 1 — bu 2 ^j 



( bu(u + l/Vb)(u-l/Vb) 
— (au — a/Vb — g'(ui)(u — l/Vb) 
—u (au — q/ Vb — g'(ui) (u — 1/ Vb) 
-b(u + l/Vb)(u-l/Vb) 

1 bu{u + l/Vb)(u-l/Vb) ^ 
(g'( Ul )-a)(u-l/Vb) 
u(g'( Ul )-a)(u-l/Vb) 
-b{u + l/Vb){u-l/Vb) 



+ o(u - Mi) 



+ o(u - Mi). 
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Dividing a common factor u—ui = u—l/Vb to these homogeneous coordinates, 
we have 

qi(w) = (bu(u + l/Vb),g'( Ul ) - a, u(^(«i) - a), -6(u + l/v / 6)) T + o(l) 
Therefore, 

Jim qi («) = (2, g'( Ul ) - a, (g'( Ul ) - a)/Vb, -2Vbf . 



It follows that qi(wi) is a finite point in 



Hence, Vi is an open curve 
in PR" 5 , since it is a continuous curve that extends to infinity in opposite 
directions with an asymptote line. 

In the third case of a < 0, it can be proved similarly that qo(«i) is a finite 
point in AIR 3 and qi(«i) is the only infinite point on V±. Therefore, in this 
case Vi is also an open curve. Finally, in all the three subcases (i.e., a = 0, 
a > and a < 0), we can show similarly that the other component V2 of the 
QSIC of A and B is also open in PR 3 . Hence, the QSIC of A and B has two 
open components in PR 3 . An example of such a QSIC is shown in Figure 7. 
This completes the proof of Theorem 6. 
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4 Classification of singular but non-planar QSIC 



4-1 [211].' /(A) = has one real double root and two other distinct roots 



Theorem 7 fl211] 3 j Given two quadrics A:X T AX = and B:X T BX = 0, 
if /(A) = has one double real root and two distinct real roots with the Segre 
characteristic [211], then the only possible index sequences of the pencil XA — B 
are (2H_2|3|2), (2» + 2|3|2), (1H_1|2|3) and (1» + 1|2|3). Furthermore, 

(1) (Case 5, Table 1) when the index sequence is (2ft_2|3|2) or (2ft + 2|3|2), 
the QSIC has one closed component with a crunode; 

(2) (Case 6, Table 1) when the index sequence is (1M_1|2|3), the QSIC has a 
closed component plus an acnode; 

(3) (Case 7, Table 1) when the index sequence is (1?? + 1|2|3), the QSIC has 
only one real point, which is an acnode. 

Proof Suppose that /(A) = has one double real root Ao and two distinct real 
roots Ai and A2 with the Segre characteristic [211]. First, it is easy to check 
that the only possible index sequences are (2??_2|3|2), (2?? + 2|3|2), (1??_1|2|3) 
and (1H + 1|2|3). 

By setting B — X$A to be B, we can transform the double root Ao into 0. 
With a further projective transform to A, we may assume that < Ai < A 2 . 
According to Theorem 1, and wlog, assuming e — 1, the two quadrics can be 
reduced simultaneously to the following forms: 



^0 10 ^ 
10 
ei 

\0 e 2 J 



B' = 




1 
eiAi 










\ 



£2X2 ) 



where £i )2 = ±1. By swapping the position a[ tl and a 44 , as well as b[ tl and 
6 4 4 , we obtain 



A' 



1 ' e 2 0^ 




1 £2X2 





o\ 


1 







1 







, B' = 








si 







eiAi 





v 1 y 




I 








There are now two cases to consider: (i) det(A') > and (ii) det(A') < 0. In 
case (i) (det(A') > 0), we have £i£ 2 = — 1- Because Id(oo)= index(A')= 2 and 
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the index jump of index function at Ao = is 0, the associated index sequence 
is (2H_2|3|2) or (2H + 2|3|2). Note that (2H_2|3|2) is equivalent to (2H + 2|1|2), 
and (2H + 2|3|2) is equivalent to (2H_2|1|2). 

By setting ^4^- A' - B' to be A', we obtain 



A' 



( r- Al-A 2 

£ 2 2 











\ 




£2X2 








(A 





-1 





r- Al+A2 









1 
















, B' = 
















A2— Ai 
2 















£1X1 





I 


r Ai+A 2 
ci 2 








/ 




I 











(12) 



Clearly, the quadric B' is a cone passing through the point (0, 0, 0, 1) T . Since 
A 2 — Ai > 0, A' is an ellipsoid if £\ = —1 and £ 2 = 1, and is a two-sheet 
hyperboloid if £\ = 1 and £ 2 = — 1. In both cases A' passes through the 
point (0,0,0, 1) T According to Eqn. (12), A' and B are in one of the two 
cases shown in Figure 8. Thus, the QSIC is a singular quartic having one 
component with a crunode. Because the QSIC is contained in the ellipsoid or 
two-sheet hyperboloid A', it is a closed curve in PR 3 . This proves the first 
item of Theorem 7. 





Fig. 8. Two cases of the QSIC having a 
crunode. 



Fig. 9. The case of the QSIC having an 
acnode. 



In case (ii) (det(A') < 0): £\ and e 2 have the same sign and Id(oo) = Id(A') = 
1 or 3. Also, the index jump at Ao = is because the size of its Jordan block 
associated with A is 2. Therefore, the associated index sequence is (1??_1|2|3) 
or (1H+1|2|3). 

If £\ = e 2 = —1, the index of XoA' — B' = —B' is 2. Thus the index sequence 
is (3?? + 3|2|l), which is equivalent to (1??_1|2|3). In this case the quadric B' 
is a cone with the |/-axis as its central axis, and the quadric A' is an elliptic 
paraboloid also with y-axis as its central axis. It is then easy to verify that 
the the QSIC has a closed component plus an acnode, as shown in Figure 9. 
This proves the second item of Theorem 7. 

If £\ = e 2 = 1, the index of XqA' — B' = —B' is 0. Thus the index sequence 
(1??1|2|3) specializes to (1H + 1|2|3). In this case, B' is a positive semi-definite; 
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thus, the QSIC of A' : X T A'X = and B' : X T B'X = has only one real 
point (0,0,0, 1), which can be verified to be an acnode. This proves the last 
item of Theorem 7. 

Theorem 8 ( / [211] 1 : Case 8, Table 1) If /(A) = has one double root and 
a pair of complex conjugate roots, then the only possible index sequence is 
(2ll_2) (or its equivalent (2ll + 2)) and in this case the QSIC comprises one 
open component with a crunode. 

Proof Suppose that the pair of complex conjugate roots are a ± bi and the 
double root is Ao- Setting B to be (B — aA)/b, we transform the roots a±bi 
to ±i. By Theorem 1, the two quadrics can be reduced to the following forms: 



A' 



'o e 0^ 




( 


eA 0^ 


e 




eA 


e 




, B' = 




1 







0-10 


v o 1 0) 




I 


lj 



(13) 



where e = ±1. Setting B' — X A to be B', and then swapping A[ v A' 44 , as 
well as B' 11 , B' 44 , A' and B' are transformed to 



A' 



10 
e 
10 
OeOO 



B' 



1 





-Ao 








e 








-Ao 





-1 


















Denote k 
with 



1 / \Jl + Ao- Applying the congruence transformation C = P T DP 



(\ A;A 0^ 



A' 



k 



/ 



10 
00 
\0 1 

simultaneously to A' and £>', we obtain the transformed A' and B' as 

/nn u n\ /} o o 0^ 

OeOO 





e 
k 2A; 2 A 
OeOO 



B' 



0-10 
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The quadric B' is a cone and therefore can be parameterized by 

r(u, v) = g(-u) + v h(w), 

where g(u) = (1-eu 2 , 2u, 1+eu 2 , 0) T and h(u) = (0,0,0, 1) T . Substituting 
r(u,v) into X T A'X = 0, we obtain a quadratic equation whose two solutions 
are v — oo, which is trivial, and t> = — c (m)/(2c!(-u)). Substituting the latter 
solution of v into r(u,v) yields the parameterization of the QSIC, 

p(«) = (2eu(l - eu 2 ), 4m 2 , 2eu(l + eu 2 ), k(l - u A ) + A; 2 A (1 + eu 2 ) 2 ) T 

(14) 

From p(u), we see that the QSIC passes through the point po = (0,0,0, 1) T 
twice, with u = or u = oo. Hence, po is a singular point of the QSIC. 
Furthermore, it is easy to verify that po is a crunode. 

In the following we will show that the QSIC has two open branches intersecting 
at the crunode. Consider the intersection of the QSIC with the plane w — 0. 
The last component w(u) of p(u) in Eqn. (14) is a quadratic polynomial in 
u 2 , whose two zeros are 

u 2 = (-ek\ + l)/(k\ -l) 

and 

u 2 = (-ek\ - l)/(k\ - 1). 

Recall that k = l/\Jl + A 2 ,, it is straightforward to verify that w(u) has two 

real zeros in u. These two real zeros are u± i2 = ±(A + \Jl + A 2 ,) when e = 1 or 
u 1;2 = ±1 when e = — 1. We observe that «i and m 2 have opposite signs and 
p(-Ui) and p(t*2) are f wo distinct points. 

Now we are going to show by contradiction that the QSIC cannot be closed. 
Assume that the QSIC is closed, i.e., there is an affine realization of PR 3 in 
which the QSIC is compact. Note that the plane w = is not necessarily 
the plane at infinity in this affine realization. Then the QSIC has a topology 
shown in Figure 10, having two closed loops joining at the crunode, i.e., like 
the figure of "8" . Since the crunode corresponds to two parameter values and 
oo of u under the parameterization p(u) in Eqn. (14), the two loops must be 
parameterized over the positive interval u G (0, +oo) and the negative interval 
u G (—00,0), respectively. Now consider again the intersection of the QSIC 
with the plane w — 0. If the plane w = intersects any loop of the QSIC, 
say the loop defined over the positive interval, there must be at least two 
intersection points, which should be given by two positive values of u through 
p(u). However, from the preceding discussions we know that there are only two 
intersections between the QSIC: p(u) and the plane w — 0, which correspond 
to one positive value and one negative value of u. This is a contradiction. 
Hence, there is no finite loop of the QSIC in any affine realization of the 
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projective space. That is, the QSIC has one open component with a crunode. 
(An example of such a QSIC is shown in Figure 11.) This completes the proof 
of Theorem 8. 



4-2 [22]: /(A) = has two double roots 

Theorem 9 ([22]2: Case 9, Table 1) If f(X) = has two real double roots 
with the Segre characteristic [22], then the only possible index sequences are 
(2ll_2ll_2) and (2ll_2ll + 2) , the QSIC comprises a real line and a space cubic 
curve intersecting at two distinct real points for both sequences. 

Theorem 10 (^22] : Case 10, Table 1) If /(A) = has two pairs of identical 
complex conjugate roots with the Segre characteristic [22], then the index se- 
quences are (2) and the QSIC comprises a real line and a space cubic curve 
that do not intersect at any real point. 

4-3 [31]: /(A) = has one real triple root and one real simple root 

Theorem 11 (Case 11, Table 1) If f(A) = has one triple root and one 
simple real root with the Segre characteristic [31], then the index sequence is 
(1M + 2|3) and the QSIC comprises a closed component with a real cusp. 

4-4 [4]-' /(A) = has one real quadruple root 

Theorem 12 (Case 12, Table 1) If /(A) = has one quadruple root with the 
Segre characteristic [4], then the index sequence is (2llll_2) or its equivalent 
form (2llU + 2), and the QSIC comprises a real line and a real space cubic curve 
tangent to each other at a real point in this case. 




Fig. 10. The hypothetical topological 
shape of a QSIC. 



Fig. 11. A QSIC having one open compo- 
nent with a crunode. 
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5 Classification of planar QSIC 

5.1 [(11)11]: /(A) = has one real double root and two other distinct roots 

Theorem 13 f[(ll)ll] 3y ) ///(A) = has one double real root and two distinct 
real roots with the Segre characteristic [(11) 11], then there are only five dif- 
ferent possible index sequences and these index sequences correspond to four 
different QSIC morphologies as follows: 

(1) (Case 13, Table 2) (2||2|1|2) - two real closed conies intersecting at two 
distinct real points; 

(2) (Case 14, Table 2) (1||3|2|3) - two real conies not intersecting at any real 
points; 

(3) (Case 15, Table 2) (1||1|2|3) - two imaginary conies intersecting at two 
distinct real points; 

(4) (Case 16, Table 2) (0||2|3|4) or (1||3|4|3) - two imaginary conies not 
intersecting at any real points. 

Theorem 14 f[(ll)ll] ly ) If /(A) = has a real double root A and a pair of 
complex conjugate roots with the Segre characteristic [(11)11], then the possible 
index sequences of the pencil \A — B are (1||3) and (2||2). Furthermore, 

(1) (Case 17, Table 2) when the index sequence is (1||3), the QSIC comprises 
of two conies, one real and one imaginary; 

(2) (Case 18, Table 2) when the index sequence is (2||2), the QSIC comprises 
of two real conies which cannot both be ellipses simultaneously in any 
affine realization o/PIR 3 . 



5.2 [(111)1] 2 ." /(A) = has one real triple root and a real simple root 

Theorem 15 ( / [(lll)l] 2y ) // f(A) = has one triple root and a simple root 
with the Segre characteristic [(111)1], then the only possible index sequences 
are (1|||2|3) and (0|||3|4). Furthermore, 

(1) (Case 19, Table 2) when the index sequence is (1|||2|3), the QSIC is a 
real conic counted twice; 

(2) (Case 20, Table 2) when the index sequence is (0|||3|4), the QSIC is an 
imaginary conic counted twice. 
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5.3 [(21) 1] 2 ■' /(A) = has one real triple root and a real simple root 

Theorem 16 ([(21)l] 2 j // f(A) = has one triple root and a simple root 
with the Segre characteristic [(21)1], then the only possible index sequences 
are (1H_|2|3) and (1» + |2|3). Furthermore, 

(1) (Case 21, Table 2) when the index sequence is (1??_|2|3), the QSIC com- 
prises two real conies tangent to each other at one real point; 

(2) (Case 22, Table 2) when the index sequence is (1?? + |2|3), the QSIC com- 
prises two imaginary conies tangent to each other at one real point. 



5.4 [2(11)] 2 .' /(A) = has two real double roots 

Theorem 17 If /(A) = has two double roots with the Segre characteristic 
[2(H)], then the only possible index sequences are (2ft_2||2) (or its equivalent 
form (2U + 2\\2)), (1H+1||3) and (1H_1||3). Furthermore, 

(1) (Case 23, Table 3) when the index sequence is (2??_2||2), the QSIC con- 
sists of a real conic and two real lines which intersect pairwise at three 
distinct real points; 

(2) (Case 24, Table 3) when the index sequence is (1^_1||3), the QSIC con- 
sists of a real conic and a pair of complex conjugate lines. The conic and 
the pair of lines do not intersect; 

(3) (Case 25, Table 3) when the index sequence is (1?? + 1||3), the QSIC con- 
sists of an imaginary conic and a pair of complex conjugate lines. The 
conic and the pair of lines do not intersect. 



5.5 [(31)]i : /(A) = has one real quadruple root 

Theorem 18 // / (A) = has one quadruple root with the Segre character- 
istic [(31)], then the only possible index sequences are (2lll_\2) and (1M + |3). 
Furthermore, 

(1) (Case 26, Table 3) when the index sequence is (2lU_\2), the QSIC consists 
of a real conic and two real lines, and these three components intersect at 
a common real point; 

(2) ( Case 27, Table 3) when the index sequence is (1M + |3), the QSIC consists 
of a real conic and a pair of complex conjugate lines, and these three 
components intersect at a common real point. 
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5.6 [(11)(11)]; /(A) = has two double roots 



Theorem 19 ([(ll)(ll)] 2y > Iff(X) = has two real double roots with the Segre 
characteristic [(11) (11)] , then the only possible index sequences are (2||2||2), 
(0||2||4) and (1||1||3). Furthermore, 

(1) (Case 28, Table 3) when the index sequence is (2||2||2) 7 the QSIC consists 
of four real lines, and these four lines form a quadrangle in PM 3 ; 

(2) (Case 29, Table 3) when the index sequence is (0||2||4) 7 the QSIC consists 
of four imaginary lines and has no real point; 

(3) (Case 30, Table 3) when the index sequence is (1||1||3) ; the QSIC consists 
of two pair of complex conjugate lines, with each pair intersecting at a 
real point. 

Theorem 20 fl(ll)(ll)] : Case 31, Table 3) If /(A) = has two identical 
pairs of complex conjugate roots with the Segre characteristic [(11) (11)], the 
only possible index sequences are (2), and in this case the QSIC comprises two 
non-intersecting real lines and two non-intersecting imaginary lines. 

5.7 [(211)]: /(A) = has one real quadruple root 

Theorem 21 If /(A) = has a quadruple root with the Segre characteris- 
tic [(211)], then the only possible index sequences are (2??_||2) and (lll_\\3). 
Furthermore, 

(1) (Case 32, Table 3) when the index sequence is (2ll_\\2) , the QSIC consists 
of a pair of intersecting real lines, counted twice; 

(2) (Case 33, Table 3) when the index sequence is (1M_||3), the QSIC consists 
of a pair of conjugate lines, counted twice. 

5.8 [(22)]: /(A) = has one real quadruple root 

Theorem 22 // /(A) = has a quadruple root with the Segre characteris- 
tic [(22)], then the only possible index sequences are (2ll_U_2) and (2U_U + 2) 
Furthermore, 

(1) (Case 34, Table 3) when the index sequence is (2ll_U_2), the QSIC con- 
sists a real double line and two other non-intersecting imaginary lines. 
The two imaginary lines do not form a complex conjugate pair. 

(2) (Case 35, Table 3) when the index sequence is (2ll_U + 2) , the QSIC con- 
sists a real double line and two other non-intersecting real lines. Each of 
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the latter two lines intersects the real double line. 



6 Classification by signature sequences 

Through the above analysis, we have put different QSIC morphologies in cor- 
respondence to different characterizing conditions, given by conditions in The- 
orem 4 through Theorem 22. Hence, we conclude that all these conditions are 
necessary and sufficient for the corresponding QSIC morphologies. We may 
then check these conditions to classify the QSIC of a given pair of quadrics in 
PM 3 . Based on these conditions one could compute the index sequence of two 
given quadrics for QSIC classification; however, this would involve the diffi- 
cult task of computing Jordan blocks. To avoid computing Jordan blocks, we 
convert all index sequences to their corresponding signature sequences. The 
advantage of using the signature sequence over using the index sequence is 
that we just need to compute the multiplicity of a real root and determine the 
signature of A A — B at the root; this is a far simpler than computing Jordan 
blocks. 

There are some cases which cannot be distinguished by using signature se- 
quences alone. We will show in this section that these cases can easily be 
distinguished by the fact that their corresponding minimal polynomials have 
different degrees. 

Not all the signature sequences of the 35 different QSIC morphologies are 
distinct: the three different QSIC morphologies with the Segre characteristics 
[lllljo, [22] and [(ll)(ll)]o (i.e., cases 4, 10 and 31) share the same index 
sequence (2), thus leading to the same signature sequence (2). Furthermore, 
two different index sequences (2\U_\2) and (2ll_ll + 2) (i.e., cases 26 and 35) are 
mapped to the same signature sequence (2((((1, 1))))2). Thus, in total, there 
are only 32 distinct signature sequences. In the following we explain how these 
cases can be distinguished. 

The signature sequence (2) can be given by the different Segre characteristics 
[lllljo, [22] and [(11)(11)] . Suppose that the signature sequence (2) has 
been detected, i.e., /(A) = has been found to have no real root. The case of 
[llll]o is distinguished from the other two cases by the fact that /(A) = has 
no multiple roots; this can be detected by whether the discriminant of /(A) 
vanishes, i.e., whether Disc(/) = ResA(/, fx) = 0. 

Then the case of [22]o and the case of [(ll)(ll)]o can be distinguished by the 
fact that they have different minimal polynomials. Suppose that the input 
quadrics are given in the real symmetric matrices A and B; and, wlog, assume 
that A is nonsingular. Since in the case of [22] or [(11)(11)] , /(A) is a squared 
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polynomial, we suppose 

/(A) = (a\ 2 + bX + c) 2 , 



whose square-free part is 

g(\) = aX 2 + b\ + c, 

where a,b,c G R and b 2 — 4ac < 0. Then, by Theorem 1 and the Cauchy- 
Cayley Theorem, the case of [(ll)(ll)]o occurs if g(X) annihilates A~ X B, i.e., 
g(A~ l B) = 0; otherwise, the case of [22] occurs. Note that g(X) can be 
obtained as the GCD of /(A) and /'(A). 

The remaining problem is that the two index sequences (2M_|2) and (2U_U + 2) 
are mapped to the same signature sequence (2((((1, 1))))2). For either of the 
two cases, /(A) = (A — a) 4 for some a G R, but the minimal polynomial for 
the case of (2ll_ll + 2) is g(X) = (A — a) 2 , while the minimal polynomial for the 
case of (2lll_\2) is h(X) = (A — a) 3 . Therefore, the case of (2ll_U + 2) occurs if 
A~ X B is annihilated by g(X), i.e., g{A~ 1 B) = 0; otherwise, the case of {2lll_\2) 
occurs. Note that g(A) = (A — a) 2 can be obtained without solving for the root 
a. 

Combining the preceding methods based on minimal polynomials with the 
methods described in Section 2.6 for exact computation of the signature se- 
quences, we have a complete algorithm for exact classification of QSIC mor- 
phologies. 

Example 1: Now we use a running example to show the procedure of using the 
signature sequence for QSIC morphology classification. Consider two quadrics 

A: 20 x 2 - 12 xy + 48 xz + 76 x + 16 y 2 

-16yz- 12y + 42 z 2 + 72 z + 58 = 0, 
B: 28x 2 + 16xy + 80xz + 56x + 2y 2 

+24yz + 20y + 56 z 2 + 72 z + U = 0. 

The equation of the eigenvalue curve C is 

m 4 + (-136A + 100) m 3 + (-1048 - 3612 A + 2904 A 2 ) u 2 
+ (-10000 A 3 + 22616 A 2 + 28416 A) u 
- 170528 A 2 + 170528 A 3 - 85264 A 4 = 0. 

Substituting u = in this polynomial yields 

-85264 A 4 + 170528 A 3 - 170528 A 2 = 0, 
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whose only real root is the double root A = 0. Substituting A = — 1 in the 
equation of C yields 

u 4 + 236 u 3 + 5468 u 2 + 4200 u - 426320, 

which has one sign change in its coefficients; therefore, by the Descartes 
rule, it has only one positive root. It follows that the signature sequence is 
(1, ((1, 1)), 3). By Theorem 14, the corresponding QSIC is the union of a real 
conic and an imaginary one, which is case 17 in Table 2. 



7 Conclusions 

To summarize, we have obtained the following result: 

Theorem 23 There are in total 35 different QSIC morphologies with non- 
degenerate pencils (see Tables 1, 2 and 3). The morphology of the QSIC of a 
pencil (A, B) is entirely classified by its signature sequence and the degree of 
its minimal polynomial, using only rational arithmetic computation. 

Besides used for determining the QSIC morphology for enhancing robust com- 
putation of QSIC in surface boundary evaluation, another application of our 
results is to derive simple algebraic conditions for interference analysis of 
quadrics. For arrangement computation, it is an interesting problem to clas- 
sify all possible partitions of M 3 that can be formed by two ellipsoids. It is 
also possible to apply the results here to derive efficient algebraic conditions 
for collision detection between various types of quadric surfaces, such as cones 
and cylinders, following the framework in (1471 ) . 

One could also use the idea developed here to study the classification of a 
pencil of conies in PR 2 , which would lead to a classification of QSIC with 
degenerate pencils in PR 3 . A more challenging problem is to use the signature 
sequence to classify the intersection of two quadrics in higher dimensions, PR 4 
say. Here the difficult issue is to deduce the geometry of the QSIC associated 
with each possible Quadric Quadric Pair Canonical Form, while it is should 
be straightforward to obtain the signature sequence of the normal form, based 
on the results presented in the present paper. 

Another direction of investigation would be the classification of the net of 
three quadrics in PK n . In this case, given three quadratic forms A, B and 
C, the question is how to use the invariants of the planar curve f(a,/3, 7) = 
det(aA + [3B + 7C) = to characterize the geometric properties of the net 
X T (aA + (3B + 7C)X = or the intersection of the quadrics X T AX = 0, 
X T BX = and X T CX = 0. 
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Proof of Theorem 9 Suppose that the two real double roots are Ao and Ai. 

By setting (B — X±A) to be B, we transform the root Ai to 0; the other root is 
still denoted by A . By Theorem 1, the two quadrics have the canonical forms 



A' = (O = 
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where e = ±1. Clearly ld(A') = Id(oo) = 2 and the index jumps at both 
roots A and Ai are 0. Therefore the only possible index sequence takes the 
form {2112R2} , covering the two nonequivalent index sequences (2ll_2ll_2) and 
(2»_2» + 2). Note that (2»_2»_2) is equivalent to (2» + 2» + 2), and (2»_2» + 2) 
is equivalent to (2ll + 2ll_2). 

Swapping a! AA and well as h' AA and b[ 1 , we obtain 
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Obviously, the QSIC contains the line x = y = 0. From the classification of 
QSIC by Segre characteristic in PC 3 , we know that the remaining component 
is a cubic curve, whose equation is found to be 

q(u) = (2e\ u/(e + u 2 ), 2X u 2 /(e + u 2 ), -u, l) T . 

It is easy to verify that the line and the cubic curve intersect at two distinct 
real points (0, 0, 0, 1) T and (0, 0, 1, 0) T . This completes the proof of Theorem 9. 

Proof of Theorem 10 Suppose that the two identical pairs of conjugate 
roots of /(A) = are a ± hi. In this case the only possible index sequence 
is (2). Setting (B — aA)/b to be B, we transform the roots a ± hi to ±i. By 
Theorem 1, the two quadrics can be transformed to the following forms, 



A' 



(q o \\ 

10 
10 
10 



B' 



1 ^ 
-1 
10 1 
0-11 i 
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Clearly, the QSIC contains the line z = w = 0, and the remaining component 
is a cubic curve parameterized by 

q( u ) = (-m 2 , u, -u(l + u 2 ), -(1 + u 2 ) ) T . 

Since the last component function — (1 + u 2 ) does not vanish for any real value 
of u, the space cubic curve has no real intersection with the line z = w = 0. 
This completes the proof of Theorem 10. 

Proof of Theorem 11 According to the discussion about the index sequences 
in Section 2.4, the index jump is 1 across the real root with a 3 x 3 Jordan 
block. Thus it is easy to see the only possible index sequence is (1M2|3). 

That the QSIC comprises a closed component with a real cusp follows from 
the classification of QSIC in PC 3 by the Segre characteristics. By complex 
conjugation it is easy to see that the cusp is real. Since the QSIC is contained 
in a projective ellipsoid in the quadric pencil (i.e., with the index being 1 or 
3), it is closed in PR 3 . This completes the proof of Theorem 11. 

Proof of Theorem 12 According to the discussion about index sequences 
in Section 2.4, the index jump is across the real root with a 4 x 4 Jordan 
block. Thus, the only possible index sequence is (2llll_2) or its equivalent form 
(2HH+2). 

That the QSIC comprises a line and a space cubic curve tangent to each 
other at a point follows from the classification of QSIC in PC 3 by the Segre 
characteristics. By complex conjugation it is easy to see that the line and the 
cubic are both real and their tangent point is also real. This completes the 
proof of Theorem 12. 

Proof of Theorem 13 Let the double zero be Ao- By setting B — XqA to B, 
we transform the double root Ao to 0. Let Ai ^ A 2 denote the other two roots. 
Wlog, we may assume < Ai < A 2 . By Theorem 1, the matrices A and B of 
two given quadrics have the following canonical forms, 

A' = diag(£i, e 2 , e 3 ,e A ), B' = diag(eiAi, e 2 A 2 , 0, ). 

Now we consider two cases: (i) det(A') > 0; and (ii) det(A') < 0. In case (i), 
the following two subcases need to be further distinguished: (i-a) £]£ 2 > and 
£364 > 0; and (i-b) E\S2 < and e 3 e 4 < 0. 

In subcase (i-a), since E\E2 > 0, £>' consists of a pair of complex conjugate 
planes; thus £>' intersects A' in two imaginary conies. Since e^e^ > 0, the index 
jump of Id(A) at A = is ±2. Since det(A') > 0, Id(-oo)=index(-A') = 0, 2, 
or 4. Therefore, all possible index sequences are (2||4|3|2), (2||0|1|2), (0||2|3|4), 
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or (4||2|1|0). Clearly, these sequences are equivalent; so we use (0||2|3|4) as 
the representative. Since there is a virtual quadric (i.e., one whose index is 
or 4), the QSIC has no real point. This completes the first part of item 4 of 
the theorem. 

In subcase (i-b), since £ 3 £ 4 < 0, the index jump of Id(A) at A = is 0. Since 
det(A') > 0, Id(— oo)= index(— A') = 0, 2, or 4. Thus, by a similar argument 
to case (i-a), the only possible index sequence is (2||2|1|2). By swapping e 3 
and £4 in the matrix A' if necessary, we may suppose that the quadric A' : 
X T A'X = is a one-sheet hyperboloid with the y-axis being its symmetric 
axis, i.e., x 2 — y 2 + z 2 — w 2 = 0. Recall that < Ai < A 2 . The two planes 
y = ±(Ai/A 2 )^ 2 x given by the quadric £>' : X T B'X = intersect A' in two 
ellipses intersecting at two real points, as shown in Figures .1 and .2. This 
completes the proof of item 1. 




Fig. .1. Cross section of a one-sheet hyper- Fig. .2. A QSIC comprising two closed 
boloid and a pair of planes. conies intersecting at two real points. 

Below we consider case (it) of det(A') < 0. Again we need to distinguish two 
subcases: (ii-a) E\Ei < and £ 3 £ 4 > 0; and (ii-b) S\E2 > and e 3 £ 4 < 0. 

In subcase (ii-a), since £364 > 0, the index jump of Id(A) at A = is ±2. 
Since det(A') < 0, Id(— 00)= index(— A') = 1 or 3. Therefore, the only possible 
index sequences are (1||3|2|3) and (1||3|4|3). Recall that E\E2 < 0. It is then 
easy to verify that the index sequence is (1||3|2|3) if S\ = —1 and £2 = 1, and 
is (1||3|4|3) if ei = 1 and e 2 = -1. 

When E\ = —1 and 62 = 1 (i.e., when the index sequence is (1||3|2|3)), by a 
projective transformation, we may transform A' to the unit sphere x 2 + y 2 + 
z 2 — w 2 = and £>' to two parallel planes z = ±(A 1 /A 2 ) 1 / 2 w. Since < Ai < A 2 , 
the QSIC consists of two real conies not intersecting each other in real points. 
This completes the proof of item 2. 

When £i = l and e 2 = —1 (i.e., when the index sequence is (1||3|4|3)), there is 
a virtual quadric in the pencil. Therefore the QSIC has no real points. Hence, 
the QSIC consists of two imaginary conies not having common real points. 
This proves the second part of item 4. 

In subcase (ii-b), using a similar argument, we know that the only possible 
index sequence is (1||1|2|3). The quadric A' is either two-sheet hyperboloid 
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with the z-axis being its centered axis (if £1, £2, and £4 have the same sign ) 
or the unit sphere centered at the origin (if £ 1; £ 2 and £3 have the same sign). 
The quadric B' comprises of a pair of imaginary conjugate planes intersecting 
in a real line — the z-axis, which intersects A' at two real points. Hence, the 
QSIC consists of two complex conjugate conies intersecting at two real points. 
This completes the proof of item 3. Hence, Theorem 13 is proved. 

Proof of Theorem 14 Wlog, by setting B — XqA to be B, we may assume 
the real double root A to be 0. Let the other two roots be a ± bi, b ^ 0. By 
Theorem 1, the two matrices A and B have the following canonical forms 



A' = diag 



V 




£1, £ 2 , II- B' = diag 



0, 0, 



V 




In the following we consider two subcases: (?) a 7^ 0; and (ii) a = 0. 
In case (i) of a ^ 0, by setting aA' — B' to be A', we obtain 



b 
-b 



A' = diag ^ £ia, £20, 
which is the quadric, 

£ l a 2 , 5 2« 2,2 2 n 

——or H — — y + z — w — 0. 


The quadric £>' consists of the following two planes 

fa±Va 2 + b 2 \ 
z = w. 



When £1 and £2 have the same sign, the index sequence is (1||3) (or its equiva- 
lent form (3| |1)). In this case, the quadric A' is either an ellipsoid or a two-sheet 
hyperboloid with two of its tangent planes being z ± w = 0. The quadric B' 
comprises two parallel planes perpendicular to z-axis. Wlog, we assume that 
b > 0. Then it is easy to verify that 

a + v^TF a- v^TF 

> 1 and - 1 < < if a > 0, 

b b 

< < 1 and < -1 if a < 0. 

b b 

It follows that one of the planes of B' intersects A' in an ellipse and the other 

plane does not intersect A' at any real point, as shown by the two cases in 

Figure .3. 



or 
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Fig. .3. The two cases in the proof of Theorem 14 where the QSIC has one real 
conic and one imaginary conic. 

When ex and e<i have opposite signs, the index sequence is (2||2). In this 
case the quadric A' is a one-sheet hyperboloid and B' comprises two planes 
intersecting A' in an ellipse and a hyperbola in the affine realization shown in 
Figure .4. Since the ellipse intersects the hyperbola at its two branches, any 
real plane in PK 3 intersects at least one of the two conies. 



Fig. .4. The case of QSIC having two real conies intersecting in two real points. 

Now consider case (ii) of a = 0. We transform the matrices A' and B' into the 
following forms, 



The quadric £>' comprises of two planes. When e 1 e 2 > 0, the only possible 
index sequence is (1||3) (or its equivalent form (3||1)), the quadric A' is a 
two-sheet hyperboloid, and the QSIC has one real and one imaginary conic, 
as shown in Figure .5. 

When E\E2 < 0, the only possible index sequence is (2||2), the quadric A' is a 
one-sheet hyperboloid, and the QSIC has two real conies intersecting in two 
real points as shown in Figure .6. Again, as in case (i) where the index sequence 
is (2||2), the two conic components of the QSIC cannot be both ellipses in any 
affine realization of PR 3 . This completes the proof of Theorem 14. 

Proof of Theorem 15 Clearly, all the roots of /(A) = are necessarily real 
in this case. Wlog, we may assume the triple root A to be zero. Then, by 
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Fig. .5. The case of the QSIC having two Fig. .6. The case of the QSIC having two 
conies, one real and one imaginary. real conies intersecting in two real points. 

Theorem 1, the canonical matrix forms of the given quadrics are 

A' = diag(ei, e 2 , e 3 , e 4 ), B' = diag(e 1 A 1 , 0, 0, ). 

Setting XiA' — B' to be A' and then E\B' to be B', we obtain 

A' = diag(0, e 2 Ai, e 3 Ai, e 4 Xi ), B' = diag(A x , 0, 0, ). 

The quadric A' is a cylinder (real or imaginary) with the x-axis being its 
central axis, while the quadric B' is a plane counted twice. Thus the QSIC is 
a conic counted twice, real or imaginary. 

When e 2 , £3 and £4 have different signs, the only possible index sequence is 
(1|||2|3) (or its equivalent form (3|||2|1)). In this case A! is real, and the QSIC 
is a real conic counted twice. When e 2 , £3 ad e 4 have the same sign, the only 
possible index sequence is (0|||3|4) (or its equivalent form (3|||0|1)). In this 
case A' is imaginary, and the QSIC is an imaginary conic counted twice. This 
completes the proof of Theorem 15. 

Proof of Theorem 16 The two roots of /(A) = are necessarily real. Let Ao 
denote the triple root and Ai denote the simple root. Wlog, we assume A = 
and Ai > 0. By Theorem 1, the canonical matrices of the two given quadric 
are 



/ 



A' 



£1 

Ex 



\ 



£2 



B' 







£1 







Since the index jump at a root of Jordan block of size 2 x 2 is 0, the index 
jump at the triple root is ±1. It is then easy to see that the only possible 
index sequence is (1??|2|3) (or its equivalents). Thus, we have e 2 = 1, £3 = 1. 

The quadric A', as 2e\ y + z 2 + w 2 = 0, passes through the origin (1, 0, 0, 0) T 
and its tangent plane at this point is y = 0. The quadric £>': E\ y 2 + Ai w 2 = 0, 
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comprises two planes intersecting at the line (y — 0,w — 0), which touches 
the quadric A' at the point (1, 0, 0, 0) T . Since the tangent plane y — of A 
at (1, 0, 0, 0) T is different from either of these two planes, the QS1C comprises 
two conies that are tangent to each other at the real point (1, 0, 0, 0) T . 

If the E\ = 1, B' consists of two imaginary planes, so the QSIC consists of two 
imaginary conies. In this case, the number of positive eigenvalues of —B' is 0, 
and the index sequence therefore becomes (1H + |2|3). If S\ = —1, £>' consists 
of two real planes, so the QSIC consists of two real conies. In this case, the 
number of positive eigenvalues of —B' is 1, and the index sequence therefore 
becomes (1??_|2|3). This completes the proof of Theorem 16. 

Proof of Theorem 17 The two double roots of /(A) = are necessarily 
real. Up to equivalence relations, the only possibilities for the index sequences 
are (2»_2||2), (1H_1||3) and (1» + 1||3). Note that (2H_2||2) is equivalent to 
(2?? + 2||2). Let Ai denote the root associated with the 2x2 Jordan block, and 
let A2 denote the other root. By setting B — \\ A to B, we may assume Ai = 0. 
Let A2 denote the other double root; wlog, assume A2 > 0. 

Then the canonical forms of the two quadrics are 

\ 



A' = diag || I , £ 2 , £3 




,.00, 

, B = diag I ! , e 2 X 2 , e 3 X 2 

0e,) J 



Due to the sign change rules (see Section 2.5), for the index sequence {2ll_2\ |2), 
we have e 2 e 3 < 0. For the index sequences (1??_1||3) and (1?? + 1||3), we have 
£2 = £3 = 1- Below we discuss these two cases: (i) e 2 e 3 < 0; and (ii) e 2 = £3 = 
1. 

In case (i) of 6263 < 0, setting e^A' — S\B' to A and setting B' to B, we 
obtain 



A = diag 



((oxA \ _ ((00, 

, 0, , B = diag I , £ 2 A 2 , e 3 \ 2 \ , 



Then the quadric A consists of two planes y = and y = 2X 2 x. The quadric 
B is a real cylinder with the its central axis. Clearly, B intersects 

y = in two real lines, denoted by i\ and £2, on the planes z — ±1. These two 
lines intersect at the point (1, 0, 0, 0) T . The quadric B also intersects the plane 
y = 2\ 2 x in a real conic, which intersects the two lines £1,2 in two distinct real 
points. This completes the proof of item 1. 

In case (ii) of e 2 = £3 = 1, Id(A') = 3, and the index jump at the root A2 is 2. 
Below we distinguish two subcases: (ii-a) S\ < 0; and (ii-b) S\ > 0. 
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In subcase (ii-a) of S\ < 0, the index sequence is (1M_1||3). The quadric B' is 
a real hyperbolic cylinder with x-axis as central axis and symmetric with the 
plane y — 0. Therefore it intersects the plane y = in two imaginary lines 
and intersects the plane y = 2\ 2 x in a real conic. The two lines intersect at 
the point (1, 0, 0, 0) T . This completes the proof of item 2. 

In subcase (ii-b) of E\ > 0, the index sequence is (1M + 1||3). The quadric B' is 
an imaginary cylinder; thus it intersects the quadric A! in an imaginary conic 
and two complex conjugate lines; these two lines intersect at the real point 
(1, 0, 0, 0) T . This completes the proof of item 3, and hence, Theorem 17. 

Proof of Theorem 18 Wlog, we may assume the quadruple root, denoted 
by Ai, to be 0. Then, by Theorem 1, the canonical form of the two quadrics 
are 



A' 







£l 



B' 







e l 
ei 







:-2) 



V e v \ 

The quadric B' comprises two planes y = and z — 0. The plane y = 
intersects the quadric A' along a real conic xz = —e2/eiw 2 , denoted as C. The 
plane z = intersects A' in two lines, defined by the intersection of z = and 
y 2 = —S2/siw 2 . If £2/^1 = —1, the two lines are real; if £ 2 /£i = 1, the two 
lines are imaginary. In both cases, the two lines intersect the conic C at the 
real point (1,0,0,0) T . 

On the other hand, when £2/^1 = — 1, we have index(A')= 2 and the index 
jump of Id(AA — B) at the root Ai = is 0. Therefore, the index sequence is 
(2M_|2) or its equivalent form (2M + |2). When £ 2 /^i — 1) we have index(A') = 
1 or 3, and the index jump of ld(XA — B) at Ai = is ±2. Therefore, the index 
sequence of the pencil is (1M|3). This completes the proof of Theorem 18. 

Proof of Theorem 19 Let the two roots be Ai and A 2 . By setting B' = 
B — X\A, we may assume Ai = 0. By Theorem 1, the canonical form of the 
two quadrics are 

A' = diag(£i, £ 2 , £3, £4), B' = diag(0, 0, £ 3 A 2 , £ 4 A 2 ). (.3) 
Setting \ 2 A' — B' to be A', we obtain 

A' = diag(£!A 2 , £ 2 A 2 , 0), B' = diag(0, 0, £ 3 A 2 , £ 4 A 2 ). (.4) 



We consider the following three cases: (i) £i£ 2 < and £ 3 £ 4 < 0; (ii) £i£ 2 > 
and £ 3 £ 4 > 0; and (in) (£i£ 2 )(£ 3 £ 4 ) < 0. 
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In case (i), the index jumps at Ai and A2 are both 0. Hence, the only possible 
index sequence is (2||2||2). In this case, each of A' and B' consists of a pair of 
real planes intersecting at a real line. Since the two real lines on A' and B' do 
not intersect, the QSIC consists of four real lines forming a quadrangle. This 
quadrangle can be obtained from a tetrahedron (defined by the four planes of 
A' and B') by removing two of the six sides; the two removed sides are the 
intersecting line of the plane pair A' and the intersecting line of the plane pair 
B' . This completes the proof of item 1. 

In case (ii), the index jumps at Ai and A2 are both ±2. Hence, the only possible 
index sequence is (0||2||4). In this case, each of A' and B' consists of a pair of 
complex conjugate planes intersecting at a real line. Since the two real lines 
on A' and B' do not intersect, the QSIC consists of four imaginary lines and 
has no real point. This completes the proof of item 2. 

In case (Hi), either E\E2 > and £364 < or E\E2 < and £364 > 0. Thus 
the only possible index sequence is (1||1||3) or its equivalent forms. At the 
same time, one of A' and B' is a pair of real planes and the other is pair of 
conjugate planes. Therefore, the QSIC consists two pairs of conjugate lines 
with each pair intersecting at a real point. This completes item 3, and hence, 
Theorem 19. 

Proof of Theorem 20 By Theorem 1, 



A' = diag 



((::) 




\ 




[10J 


J 



, B' = diag(-l, 1, -1, 1). 



Then A' is a hyperbolic paraboloid xy + zw = and B' is a hyperboloid 
— x 2 + y 2 — z 2 + w 2 = 0. It is easy to verify that the QSIC consists of the two 
non-intersecting real lines defined by (x + w = 0, y — z = 0) and (x — w = 0, 
y + z = 0), and two non-intersecting imaginary lines defined (x — iz = 0, 
y — iw = 0) and (x + iz = 0, y + iw = 0). In this case, since /(A) = has no 
real root, the index sequence is (2). This completes the proof of Theorem 20. 

Proof of Theorem 21 The quadruple root, denoted by Ai, of /(A) = is 
necessarily real. Wlog, we may assume Ai = 0. By Theorem 1, the canonical 
form of the two given quadrics are 



A' = diag \,e 2 ,e 3 \,B , = diag(0, e u 0, 




The quadric B' is a pair of the identical real plane y — 0. Substituting y 2 = 
in the the quadric A', we find that the QSIC is the intersection between 
y 2 = and the pair of planes z 2 = —e 3 /e 2 w 2 . Clearly, the QSIC comprises 
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two intersecting real lines, counted twice, if £2^3 = — 1, or a pair of conjugate 
lines, counted twice, if £ 2 £3 = 1- 



Clearly, the index jump of ld(XA' — B') at Ai = is or ±2. When e 3 /e 2 = — 1, 
det(A) = —6263 > 0, therefore index(A') is 0, 2 or 4. It follows that the 
index sequence is (2^_||2) or its equivalent form (2H + ||2). When £3/62 — 1, 
det(A) = —£2^3 < 0, therefore index (A') is 1 or 3. It follows that the index 
sequence is (1M_||3) or its equivalent form (1H + ||3). This completes the proof 
of Theorem 21. 

Proof of Theorem 22 The quadruple root of /(A) = is necessarily real, 
and may be assume to be 0. Since the index jump at a real root with a 2 x 2 
Jordan block is 0, the index sequence is of the form (211112). By Theorem 1, 
the canonical form of the two quadrics are 



and £>' is a pair of planes 

e x y 2 + e 2 w 2 = 0. 

When E\E2 = 1, the index sequence is (2ll_U_2) or its equivalent form (2ll + ll + 2). 
In this case B' comprises two conjugate planes y+iw = and y—iw = 0, which 
intersects A' in the real double line (y — 0, w — 0) and two non-intersecting 
imaginary lines (x — iz = 0, y — iw = 0) and (x + iz = 0, y + iw = 0). 

When e±e2 = —1, the index sequence is (2ll_U + 2) or its equivalent form 
(2ll + ll_2). In this case £>' comprises two real planes y — w = and y + w = 0, 
which intersects A' in the real double line (y — 0, w — 0) and two non- 
intersecting real lines (x — z = 0, y — w = 0) and (x + z = 0, y + w = 0). 

This completes the proof of Theorem 22. 




The quadric A' 



is the hyperbolic paraboloid 



E\xy + e 2 zw = 0, 
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